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Abstract 

This  thesis  presents  a  constrained  convex  optimization  control  design  procedure  that 
takes  advantage  of  a  decoupled  hierarchical  system  structure.  To  achieve  a  hierarchi¬ 
cal  structure,  a  system  is  decoupled  into  primary  channels  that  have  its  own  dynamics, 
control  actuator,  and  a  resulting  controller.  These  MIMO  (multiple  input/multiple 
output)  channels  are  then  decoupled  into  a  hierarchical  structure  of  SISO  (single  in¬ 
put/single  output)  loops.  The  constrained  optimization  procedure  is  then  applied  to 
each  SISO  loop  to  determine  a  solution  controller. 

The  procedure  minimizes  the  closed  loop  7^2  norm  as  the  primary  objective,  while 
presenting  the  Q-minimization  objective  as  a  secondary  option.  The  procedure  applies 
closed  loop  frequency  and  transient  response  constraints  to  the  constrained  optimiza¬ 
tion,  which  utilizes  Q-parameterization  of  the  objective  function.  The  Q-parameter 
optimization  decision  variables  are  the  coefficients  of  a  stable  orthogonal  basis.  This 
basis  is  constructed  using  a  presented  Basis  Function  Algorithm,  designed  as  a  con¬ 
sistent  procedure  for  basis  selection.  The  algorithm  is  designed  to  create  the  most 
efficient  basis  from  three  forms  of  basis  functions:  Finite  Impulse  Response  (FIR) 
functions,  Laguerre  functions,  and  Legendre  functions.  User  defined  Fixed  Pole  Model 
basis  functions  are  also  addressed. 

The  design  method  is  applied  using  the  Draper  Laboratory’s  Structural  Control 
Toolbox  (SCTB),  a  Computer  Aided  Design  (CAD)  tool.  The  design  procedure  is 
used  on  the  Draper  Small  Autonomous  Aerial  Vehicle  (DSAAV)  helicopter  as  a  design 
example.  The  solution  procedure  was  able  to  make  improvements  in  the  design  of 
the  DSAAV  forward  motion  channel.  This  design  utilized  the  modularity  of  the 
hierarchical  structure,  using  both  constrained  optimization  and  classical  controller 
designs. 

Thesis  Supervisor:  Dr.  Marc  W.  McConley 
Title:  Senior  Member  of  the  Technical  Staff 

Thesis  Advisor:  Professor  Eric  Feron 

Title:  Assistant  Professor  of  Aeronautics  and  Astronautics 
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Chapter  1 


Introduction 

1.1  Historical  Overview 

Control  engineering  is  rooted  in  two  main  methods  of  solving  problems:  classical 
methods  and  model-based  methods.  Classical  control  methods,  such  as  use  of  root- 
locus  techniques,  Bode  plots,  and  Proportional  Integral  Derivative  (PID)  controllers 
infer  closed-loop  performance  and  are  very  effective  for  single  input/single  output 
(SISO)  systems.  However,  applying  classical  control  to  multiple  input/multiple  out¬ 
put  (MIMO)  systems  is  of  limited  value  because  of  the  coupling  of  the  inputs  and 
outputs  within  the  dynamics.  Also,  classical  methods  deal  with  design  specifications 
indirectly,  requiring  the  control  engineer  to  iterate  by  trial  and  error  until  a  solution 
is  found. 

Model-based  methods,  such  as  the  “Hoo  aJid  Linear  Quadratic  Gaussian  (LQG  or 
7^2)  methods,  are  based  on  state  space  representations  of  the  system.  These  methods 
map  the  closed  loop  system  in  terms  of  a  metric  or  norm.  Often  frequency  weights 
are  added  to  the  input  and  output  signals  to  shape  the  closed  loop  response.  These 
weights  become  the  design  variables  and  are  adjusted  until  the  desired  closed  loop 
performance  is  achieved.  Although  model-based  methods  are  effective  in  controlling 
MIMO  systems,  the  adjustment  of  these  design  parameters  is  still  a  difficult  and 
indirect  way  to  design  controllers. 

With  computer  technology  becoming  faster  and  cheaper,  convex  optimization  has 
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been  researched  as  a  more  direct  method  of  controller  design.  Convex  optimization 
is  distinct  from  non-convex  optimization  because  convex  optimization  does  not  have 
local  minima.  Therefore,  any  minimum  in  a  convex  problem  is  a  global  minimum. 
Important  convex  optimization  programs  are  the  linear  program  and  the  quadratic 
program.  A  linear  program  has  a  linear  objective,  such  as  c^x,  and  linear  constraints. 
A  popular  technique  for  solving  linear  programs  is  the  Simplex  Method  [17].  A 
quadratic  program  has  a  quadratic  objective  function,  such  as  x’^Cx.  If  the  constant 
matrix  C  is  positive  definite  and  the  constraints  are  linear,  then  the  quadratic  program 
is  convex.  In  the  1960s,  Fegley  and  colleagues  [2],  investigated  the  application  of  linear 
and  quadratic  programming  to  optimal  control  problems. 

In  1976,  Youla  first  recognized  Q  (or  Youla)  parameterization  [18].  In  the  1980s, 
^-parameterization  methods  of  controller  design  were  developed,  applying  the  Q 
parameter  to  the  closed  loop.  Gustafson  and  Desoer  generalized  the  parameterization 
of  Q,  representing  Q  in  terms  of  zero  and  pole  locations  [3,  4].  Boyd  [1]  and  Polak  [15] 
used  ^-parameterization  in  a  convex  optimization  control  design.  This  was  done 
through  a  finite  approximation  of  Q  via  basis  functions.  Optimization  over  all  closed 
loops  was  reduced  to  an  optimization  over  all  stable  Q.  This  is  a  foundation  of 
constrained  convex  optimization. 

Most  recently,  constrained  optimization  was  successfully  applied  by  McGovern, 
using  'H2  and  £i  objectives  on  an  actual  hardware  system,  the  structural  control  of  a 
telescope  [9,  10].  Also,  Lintereur  has  successfully  applied  constrained  optimization, 
using  the  "Hz  objective,  on  a  spacecraft  attitude  control  system  [6,  7]. 

1.2  Contributions 

There  are  two  main  goals  of  this  thesis.  One  is  to  uncover  some  benefits  of  decoupling 
a  system  into  a  hierarchical  structure  and  applying  the  constrained  optimization  one 
element  at  a  time.  These  benefits  could  range  from  improving  optimization  efficiency 
to  taking  advantage  of  modular  control  design,  where  controllers  for  separate  modes 
(within  the  same  channel)  could  be  designed  with  completely  different  techniques. 
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The  other  goal  is  to  apply  constrained  optimization  to  an  actual  control  problem. 
The  actual  control  problem  presented  in  this  thesis  is  the  Draper  Small  Autonomous 
Aerial  Vehicle  (DSAAV),  a  helicopter. 

In  addition  to  the  above  goals,  this  thesis  makes  several  further  contributions.  One 
such  contribution  is  a  procedure,  based  on  constrained  optimization,  that  solves  for 
controllers  within  the  hierarchical  structure.  Another  contribution  is  the  application 
of  the  Structural  Control  Toolbox  (SCTB)  to  the  solution  procedure.  The  SCTB, 
developed  at  the  Draper  Laboratory,  is  a  Computer  Aided  Design  (CAD)  toolbox 
which  runs  under  MATLAB  and  can  be  very  useful  for  many  forms  of  controller 
design.  This  thesis  explains  how  to  apply  the  powerful  SCTB,  detailed  in  Chapter  5, 
and  the  benefits  of  this  relatively  new  toolbox.  Another  useful  aspect  is  an  algorithm 
for  choosing  the  most  efficient  bctsis,  within  the  solution  procedure,  based  on  the 
forms  of  basis  functions  available  in  the  SCTB.  This  algorithm  aims  at  reducing  the 
number  of  basis  functions  necessary  to  estimate  the  objective  function  in  order  to 
improve  the  speed  of  the  optimization  and  produce  lower  order  controller  solutions. 


1.3  Organization 

The  organization  of  this  thesis  falls  into  seven  chapters.  The  second  chapter  provides 
developmental  information  for  constrained  optimization.  This  starts  with  a  detailed 
explanation  of  how  to  formulate  Q-parameterization  into  a  constrained  optimization 
problem.  It  then  explains  the  different  options  for  the  portions  of  the  constrained 
optimization  problem:  choosing  an  objective  function,  applying  constraints,  and  se¬ 
lecting  a  ba;sis  function  representation  of  Q.  The  chapter  ends  detailing  the  model 
reduction  techniques  necessary  for  applying  constrained  optimization  to  higher  order 
systems. 

Chapter  3  details  the  decoupled  hierarchical  system  architecture  and  how  to  break 
a  coupled  system  into  the  hierarchical  structure.  Chapter  4  details  the  solution  pro¬ 
cedure,  complete  with  algorithms  for  choosing  objective  functions  and  setting  con¬ 
straints.  There  is  also  a  detailed  Basis  Function  Algorithm  in  this  chapter,  designed 
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to  find  the  most  efficient  basis  that  will  provide  a  feasible  solution.  Chapter  5  ex¬ 
plains  the  SCTB  and  how  to  apply  it  to  the  solution  procedure.  Chapter  6  applies 
the  solution  procedure  to  the  DSAAV  helicopter  in  a  control  design  example.  Finally, 
the  conclusions  and  suggestions  for  future  work  are  presented  in  Chapter  7. 
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Chapter  2 


Constrained  Optimization 


This  chapter  presents  a  brief  overview  of  some  concepts  necessary  for  the  development 
of  the  Constrained  Optimization  solution  procedure  in  this  thesis.  It  will  be  useful 
for  understanding  concepts  and  notation  used  in  later  chapters.  The  reader  may  skip 
this  chapter  if  these  axe  familiar  concepts  and  should  refer  back  as  necessary. 

A  general  control  optimization  problem  is  presented  based  on  the  generic  system 
in  Figure  2.1.  The  closed  loop  system  H  is  represented  by  n„  X  Uy  compensator  K 
and  riz  X  nt„  system  H  =  K).  Here,  J'liP,  K)  denotes  the  linear  fractional 

transformation  of  the  plant  P  with  the  compensator  K.  The  objective  function  that 
determines  how  the  closed  loop  performance  is  measured  is  ^{H).  The  closed  loop 
system  H  is  constrained  to  the  closed  loop  constraint  set  JCu-  Enough  information  is 
provided  to  define  the  following  general  control  optimization  problem. 

Problem  1  Given  the  closed  loop  constraint  set  ICh)  choose  the  stabilizing  compen¬ 
sator  K  that  solves  the  following  optimization  problem,  minimizing  the  closed  loop 
system  H: 

min 

Kstgbilizing 

subject  to  H  —  iF({PyK)  G  ICh 

Problem  1  is  a  non-convex  optimization  problem,  which  is  difficult  to  solve.  Q- 
parameterization  (or  Youla  Parameterization)  allows  for  parameterization  of  all  sta- 
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Figure  2.1:  General  Control  Problem 

bilizing  controllers  by  the  parameter  Q.  Q  lies  in  a  convex  set  where  H  is  linear  in 
Q,  maldng  for  a  convex  optimization  problem  that  is  easier  to  solve. 

This  thesis  will  define  constrained  optimization  as  the  important  application  of 
the  powerful  ^-parameterization  allowing  for  convex  optimization  formulation.  This 
formulation  allows  ICh  to  be  mapped  to  convex  constraints  on  Q.  Every  stabilizing 
compensator  is  related  to  Q  through  a  bilinear  map.  The  mathematical  formulation 
is  the  following  optimization  problem  [1,  7]. 

Problem  2  Given  the  closed  loop  constraint  set  K-u  and  the  set  of  Q  that  closes  loops 
in  Kh)  ICq,  choose  the  infinite  dimensional  design  variable  Q  to  solve  the  following 
optimization  problem  that  minimizes  the  closed  loop  system  H: 

min  ^(H) 

QeKg 

subject  to  Kq  =  {Q  stable  |  Ti  -1-  T2QT3  G  ICh} 

H  =  Ti+T2QT3 

The  equation  that  defines  the  closed  loop  system  H  is  the  ^-parameterization  of  H, 
later  defined  as  Equation  2.1. 

This  chapter  starts,  in  Section  2.1,  with  a  discussion  of  ^-parameterization,  the 
parameterization  technique  that  makes  constrained  optimization  possible.  The  fol¬ 
lowing  section,  Section  2.2,  defines  the  possible  objective  functions  used  to  determine 
the  performance  of  the  optimization.  Section  2.3  details  basis  functions  and  how  they 
represent  Q.  This  section  explains  each  of  the  four  basis  functions  available  to  create 
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Figure  2.2:  Q-Paxameterization. 

the  orthogonal  basis  for  the  decision  variables.  Section  2.4  describes  how  constraints 
on  the  frequency  and  time  response  are  formulated  into  constrained  optimization 
constraints.  The  final  section  in  the  chapter  explains  the  model  reduction  techniques 
necessary  for  applying  constrained  optimization  to  higher  order  systems. 


2.1  Q-Parameterization 

The  first  step  in  solving  the  constrained  optimization  problem  is  to  find  the  Q- 
Parameterization,  or  the  Youla  Parameterization  [8,  18],  of  the  closed  loop  system. 
The  closed  loop  system  is  represented  by  n«  x  Uy  compensator  K  =  J-t{KaiQ) 

Uz  X  Uto  system  H  =  .F/(P,  K)  in  Figme  2.2. 

2.1.1  Q  Representation  of  the  Closed  Loop  System 

Q-parameterization  is  achieved  by  developing  the  stabilizing  augmented  observer- 
based  controller  Ka,  with  augmented  control  input  v  (added  to  the  actuator  command 
signal)  and  augmented  observer  error  output  e.  The  new  output  e  becomes  the  input 
and  the  new  input  v  becomes  the  output  for  the  unknown  stable  and  realizable  Q 
_  Qe),  por  any  stable  Q,  this  does  not  change  the  stability  of  the  closed  loop 
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system.  Furthermore,  it  is  well  known  that  the  closed  loop  map  from  n  to  e  is  zero. 
This  allows  Q  to  appear  linearly  in  the  closed  loop  system  H: 

H  =  Ti  +  T2QT3  (2.1) 

where  Ti  is  the  open  loop  map  from  w  to  z,  T2  is  the  map  from  w  to  e,  and  T3  is  the 
map  from  v  to  z.  H  can  now  be  represented  by  Equation  2.1  as  opposed  to  the  linear 
fractional  representation  The  parameterization  in  Equation  2.1,  affine  in 

Q,  represents  all  stabilizing  controllers,  allowing  a  search  over  all  realizable  closed 
loop  systems  to  be  replaced  by  a  search  over  all  stable  Q.  Similar  detailed  discussion 
of  ^-parameterization  can  be  found  in  [1,  5,  6,  7,  8,  9,  18]. 

2.1.2  Q-Parameterization  in  State  Space 

The  plant  P  maps  the  exogenous  input  w  and  control  u  to  the  regulated  output  z 
and  measurement  y  as  follows: 


For  the  compensator  K,  the  closed  loop  system  H  is: 

H  =  Pn  +  PuK{I  -  P22K)-^Pn  =  W  K)  (2.3) 

For  ^iiT  to  be  internally  stabilizing,  H  can  be  replaced  by  the  ^-parameterization  of 
Equation  2.1. 

The  Q  parameter  can  be  derived  from  any  nominally  stabilizing  controller.  The 
state  space  description  of  the  discrete  plant  P  is: 

x[k  -1-1]  A  Bi  B2  a:[fc] 

z[k]  =  Cl  Dll  Di2  «;[A:]  •  (2.4) 

y[k]  C2  D21  D22 

A  controller  gain  matrix  F  and  observer  gain  matrix  G  are  chosen  to  place  the  eigen¬ 
values  of  {A  —  B2F)  and  {A  —  GC2)  inside  the  unit  circle  to  create  a  stabilizing 
model-based  controller.  To  augment  the  controller  with  Q  as  in  Figure  2.2,  an  input 
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V  is  added  to  the  actuator  command  signal.  The  output  measurement  residual  e  and 

V  are  described  in  the  following  equations: 

u[A:]  =  Qe[k] 

u[A:]  =  -Fx[k]  +  u[fc]  (2.5) 

e[k]  =  y[k]  —  C2x[k]  —  D22u[k] 

where  x  represents  the  state  estimate.  A  stable  Q  will  not  affect  the  stability  of 
the  closed  loop  system.  After  including  the  contribution  of  v  to  the  state  estimate 
equation,  the  state  space  representation  of  the  augmented  controller  Ka  is: 


x[k  + 1] 
«[fc] 
e[k] 


A-B2F-GC2  +  GD22F  G  B2-GD2 


—C2  ~  D22F 


where  the  space  of  cill  stabilizing  controllers  is  spanned  hy  K  =  if  Q  is 

stable. 

Algebraically  eliminating  u  and  y  in  a  lower  linear  fractional  transformation  of  P 
and  Ka  results  in  the  following  state  space  representation  of  T: 


x[k  + 1] 
x[k  + 1] 
z[k] 
e[k] 


A  +  B2F 


A+GC2  B1  +  GD21  0 


Cl  +  D12F  —D12F 


where  the  state  estimate  error  x  =  x  —  x.  The  state  estimate  error  is  uncontrollable 
from  V  and  the  state  x  is  unobseirvable  from  e.  Therefore,  the  transfer  function  from 
t;  to  c  is  zero.  Figure  2.3  shows  the  relationship  between  T  and  Q.  The  closed  loop 


from  w  to  z  is: 


H{z)  =  Ti  +  T2Q{I  -  =  Te{T,  Q) 


where  T  =  .  Since  T4,  the  transfer  function  from  v  to  e,  is  zero,  Equa- 

[  Ta  T4  J 

tion  2.8  reduces  to  Equation  2.1,  where  the  elements  of  T  are: 


A-^  B2F  — 

0  H"  GC2 


J5  1 

,  '  ,  C1  +  D12F  -DiiF  ,  Dn 

Bi  +  GD21 
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Figure  2.3:  T  and  Q  Depiction  of  System 

Tz  =  (  A  +  B2F,  B2,  Cl  +  D12F,  Di2  )  (2.10) 

T3  =  {  A  +  GC2,  Bi  +  GD21,  C2,  D21  )•  (2.11) 

If  the  open  loop  plant  is  stable,  the  controller  gain  matrix  F  and  observer  gain  matrix 
G  can  go  to  zero.  This  results  in  Ti  =  Pn,  T2  =  P121  «md  T3  =  P21. 

2.1.3  Infinite  Dimensional  Representation  of  Q{z) 

Q  can  be  represented  by  an  infinite  number  of  basis  functions  and  coefficients.  This 
infinite  representation  of  Q  is  as  follows: 

Q{^)  =  flF„{z)<l>„  (2.12) 

n=0 

where  Fn{z)  is  a  complete  set  of  basis  functions  and  <f)n  represents  the  free  basis 
function  coefficients.  To  be  complete,  Fn{z)  must  span  the  space  of  all  stable  real¬ 
izable  transfer  functions.  The  relationship  between  this  representation  of  Q,  basis 
functions,  and  optimization,  including  the  finite  representation  of  Q{z),  is  discussed 
in  Section  2.3. 


2.2  Objective  Functions 

The  objective  function  defines  how  the  constrained  optimization  measures  system 
performance.  The  choice  of  objective  function  can  radically  alter  how  the  problem 
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is  formulated  and  the  success  of  the  resulting  controller.  There  are  various  types  of 
objective  functions,  and  the  SCTB  looks  at  three:  'H.2  minimization,  Q-minimization, 
and  minimization.  This  thesis  will  investigate  the  first  two:  7^2  minimization  and 
^-minimization. 

2.2.1  ^.2  Minimization 

The  most  effective  way  to  initially  design  a  constrained  optimization  controller  is 
through  use  of  an  'H2  objective  function  [5].  This  objective  function  will  minimize 
the  7^2  norm  of  the  transfer  matrix  from  disturbance  and  noise  inputs  to  error  and 
control  outputs.  Looking  back  to  Figure  2.2,  there  are  two  kinds  of  plant  inputs:  the 
exogenous  input  iw,  containing  disturbances  and  noises,  and  the  controller  output  u. 
Likewise  there  are  two  kinds  of  plant  outputs:  the  regulated  output  z,  containing 
errors  and  control  signals,  and  the  controller  input  y. 

Q-Parameterization  of  'H2  Norm 

One  of  the  most  popular  standards  of  measurement  of  a  model-based  system  is  the 
system  7^2  norm.  This  norm  is  the  root  mean  squared  (RMS)  value  of  the  system 
output  to  a  given  unit  variance  (Gaussian)  white  noise  input.  The  H2  norm  of  a 
MIMO  system  is: 

n*  nw  00 

«=i  i=i  fc=o 

where  h[k]  represents  the  closed  loop  impulse  response. 

The  Q-parameterization  representation  of  the  7^2  norm  is  based  on  the  following 
Q-parameterization  representation  of  the  closed  loop  impulse  response: 

flu  JV— 1 

hij[k]  =  fn4>pqn  *  ^3,gi)[fc]  (2*14) 

p=l  g=l  n=0 

where  t2,iji  and  tz,ij  are  the  impulse  responses  of  the  (i,  j)  elements  of  Ti,  r2, 
and  Tz  respectfully  and  the  fn4>pqn  term  provides  the  impulse  response  of  Q  for  iV 
finite  number  of  ba.sis  functions  /„  and  free  coefficients  <f>pqn’ 
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■^2  Objective  Function 


The  7^2  suboptimal  objective  function  is  the  following  quadratic 

vector  representation  of 

function  of  <l>pq,  the 

riu 

min  £  £  4>lqMpq(t)pq  + 

^  P=1  ^=1 

(2.15) 

where  the  terms  Qpq^  and  Mpq  are  defined  as  follows: 

*  [fo  fl  ’  *  *  /iVi  ] 

(2,16) 

nz  nw 

9pq  ~ 

i=l  j=l 

(2.17) 

nz  nw 

t=l  J=1 

(2.18) 

Once  constraints  are  added,  the  7^2  minimization  constrained  optimization  problem 
is  represented  by  the  convex  quadratic  program  [5]: 


min  g^<f> 

4> 

subject  to  A(j)  <  6 

where  A  and  b  are  derived  from  user  defined  constraints.  Constraints  axe  formulated 
in  Section  2.4.  Similar  explanations  of  the  7^2  objective  function  are  found  in  [5,  7,  9]. 

2.2.2  Q-Minimization 

The  Q-minimization  technique  attempts  to  improve  upon  an  existing  stabilizing  nom- 
inal  controller  K„om-  This  form  of  constrained  optimization  designs  an  augmentation 
to  the  nominal  controller.  This  is  illustrated  in  Figure  2.4,  where  the  nominal  plant 
P  and  Knom  are  closed  to  form  the  nominal  closed  loop  Hnom‘  The  constrained 
optimization  design  is  the  augmentation  controller  Kaug- 

With  the  stability  of  the  nominal  closed  loop  system  ffnom  as  a  given,  an  observer 
based  Q-paxameterization  can  occur  without  the  observer  and  state  feedback  gains 
previously  required  to  stabilize  the  system.  This  stabilizing  controller  Kg  is  then 
combined  with  Q  to  complete  the  augmentation  controller. 

=  (2.19) 
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Figure  2.4:  Augmented  Controller  Design 

The  nominal  controller  is  then  added  to  complete  the  iinal  augmented  controller  K. 

K  —  Knom  +  Kaug  (2.20) 

It  cannot  be  stressed  enough  that  a  good  K„om  is  necessary.  As  shown  in  Equa¬ 
tion  2.20,  any  undesirable  characteristics  of  Knom  will  very  likely  be  retained  in  the 
final  design. 


(^-Minimization  Objective  Function 


A  good  Knom  is  required  because  Q-Minimization  tries  to  minimize  the  norm  of 
Q,  as  opposed  to 7^2  minimization,  which  attempts  to  minimize  the  7^2  norm  of  the 
closed  loop  system.  This  way,  the  modifications  to  Knom  are  penalized,  with  the  worst 
case  scenario  of  Q  =  0  resulting  in  an  unaugmented  Knom  as  the  fin^  controller.  For 
orthonormal  basis  fimctions,  the  7^2  norm  of  Q  is  defined  as  follows: 


I|Q(^)I|2 


riu 


*=1  j=l  n=0 


(2.21) 


Now,  the  suboptimal  Q-minimization  objective  function  is: 


nun^^^ 


(2.22) 
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with  <f)  being  the  vector  representation  of  <^,yn*  Similar  explanations  of  the  Q-minimi- 
zation  objective  function  axe  found  in  [5,  7,  9]. 


2.3  Basis  Functions 

In  the  Q-paxameterization-based  optimal  control  formulation  of  Problem  2,  Q  repre¬ 
sented  the  infinite  dimension  design  variable.  As  depicted  in  Equation  2.12,  Q  can  be 
defined  as  a  function  of  coefficients  (j)  and  basis  functions  Fn{z).  This  can  be  easily 
approximated  to  a  finite  representation  of  Q  with  a  finite  number  of  basis  functions. 

There  axe  several  different  forms  of  basis  functions.  This  thesis  will  focus  on 
the  four  kinds  of  basis  functions  that  axe  available  in  the  SCTB:  Finite  Impulse 
Response  (FIR),  Laguerre  functions,  Legendre  functions,  and  the  user  defined  Fixed 
Pole  Model  (FPM)  [5].  Further  discussions  on  FIR  and  Lagaurre  basis  functions  axe 
found  in  [5,  7,  9] 

2.3.1  Finite  Representation  of  Q{z) 

To  make  the  optimization  problem  tractable,  it  is  necessaxy  to  make  a  finite  approx¬ 
imation  of  Q.  A  finite  approximation  of  Equation  2.12  is  defined  below: 

Q{z)  Fr,{z)<f>„  (2.23) 

n=0 

where  Fn{z)  is  the  sequence  of  N  stable  basis  functions.  Now,  the  corresponding 
set  of  free  coefficients,  becomes  the  decision  variables  for  the  objective  functions 
described  in  Section  2.2.  It  should  be  noted  that  if  F„(z)  is  orthonormal  and  N  goes 
to  infinity.  Equation  2.23  can  represent  any  stable  transfer  function. 

2.3.2  Finite  Impulse  Response 

Simplicity  makes  the  FIR  the  most  popular  basis.  The  FIR  formulation  assumes  the 
impulse  response  of  Q  is  zero  after  N  time  steps.  The  FIR  representation  of  Q{z)  is: 

Q(z)  =  E  (2.24) 

n=0  ^ 


30 


This  basis  is  orthonormal,  satisfying  the  following  relationship: 

r 

oo  1  for  %  —  7 

/.  ■  (2-25) 

k=o  0  for  i  ^  j 

where  the  time  domain  form  of  the  functions  are  /,[fc]  =  5{i  —  k). 

The  FIR  basis  is  best  used  for  modeling  fast,  highly  damped  systems.  However,  if 
the  optimal  Q  contains  lightly  damped  modes  or  very  slow  dynamics,  the  FIR  basis  is 
less  effective.  Slow  systems  can  require  an  inefficiently  high  number  of  basis  functions 
iV,  limiting  the  effectiveness  of  the  FIR  basis.  It  is  suggested  that  FIR  basis  functions 
are  most  effective  in  combination  with  other  forms  of  basis  functions  [5]. 


2.3.3  Laguerre  Functions 


Laguerre  functions  axe  a  generalization  of  the  FIR  basis.  The  difference  between  the 
two  basis  forms  is  Laguerre  functions  have  their  poles  placed  at  a  specified  time  scale 
a,  given  that  |a|  <  1.  The  FIR  has  all  poles  placed  at  the  origin. 

The  discrete  Laguerre  basis  representation  of  Q{z)  is: 


^n- 


(2.26) 


Notice  that  for  a  =  0,  the  Laguerre  basis  is  equivalent  to  the  FIR  basis  of  Equa¬ 
tion  2.24.  Every  transfer  function  has  an  associated  optimal  time  scale  a.  Similar  to 
FIR,  the  Laguerre  basis  functions  are  orthonormal  to  each  other  and  restrict  all  of 
the  poles  to  one  location,  but  the  user  specifies  the  location  through  the  choice  of  a. 


Selection  of  the  Time  Scale  a 

The  number  of  basis  functions  N  necessary  to  adequately  approximate  the  optimal  Q 
is  directly  related  to  the  choice  of  a.  A  poorly  chosen  a  will  result  in  a  very  large  basis, 
where  a  well  chosen  a  will  keep  the  basis  size  small.  A  well  chosen  a  will  place  it  close 
to  the  dominant  dynamics  of  the  optimal  Q,  However,  selection  of  the  best  time  scale 
and  corresponding  number  of  basis  functions  is  hardly  intuitive  and  often  requires 
several  design  iterations.  The  SCTB  suggests  choosing  a  time  scale  with  a  frequency 
in  a  region  where  the  problem  is  highly  constrained  in  the  frequency  domain  [5]. 
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2.3.4  Legendre  Functions 


Legendre  functions  axe  an  orthonormal  set  of  basis  functions  with  poles  spread  over 
a  range  of  frequencies.  This  differs  from  the  FIR  and  Laguerre  basis,  which  have 
their  functions  located  at  the  same  pole.  The  discrete  Legendre  basis  representation 
of  Q{z)  is: 


Q{z) 


(2.27) 


where  Cjt  =  ae  (*+o-®)  and  |al  <  1. 


Selection  of  the  Time  Scale  a 

As  with  Laguerre  functions,  the  time  scales  o,-  determine  the  location  of  the  basis 
poles.  Poles  of  Legendre  functions  are  distributed  across  the  frequency  band  at  0.5, 
1.5,  2.5,  3.5,  . . .,  (iV  —  0.5)  times  the  frequency  associated  with  the  chosen  time  scale 
a.  For  example,  if  the  chosen  a  has  a  frequency  of  10  Hz,  the  corresponding  Legendre 
basis  would  have  poles  located  at  5  Hz,  15  Hz,  25  Hz,  35  Hz,  . . .,  (ION  —  5)  Hz. 
The  SCTB  User’s  Guide  suggests  choosing  a  time  scale  frequency  no  higher  than  the 
frequency  of  the  lowest  expected  (non-integrator)  controller  pole  or  zero  [5].  With 
their  dynamics  distributed  across  a  range  of  frequencies,  Legendre  functions  provide 
for  a  very  flexible  basis.  This  flexibility  makes  the  Legendre  functions  a  great  initial 
guess  for  a  basis. 

2.3.5  Fixed  Pole  Model 

The  fixed  pole  model  (FPM)  basis  allows  the  user  to  determine  the  pole  locations  of 
the  basis.  The  SCTB  gives  the  user  the  opportunity  to  define  FPM  poles  based  on  the 
plant  poles  and  the  Q  poles,  from  the  most  recent  design  of  Q  [5].  More  information 
on  how  the  SCTB  allows  the  user  to  use  the  plant  poles  and  Q  poles  can  be  found  in 
Chapter  5  or  [5]. 
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2.4  Closed  Loop  Constraints 


After  defining  the  objective  function  and  basis,  the  remaining  task  t'o  complete  the 
formulation  of  the  constrained  optimization  problem  is  to  define  constraints.  For  the 
purposes  of  this  thesis,  constraints  are  either  applied  to  the  closed  loop  frequency 
response  or  the  closed  loop  time  response. 

2.4.1  Frequency  Response  Constraints 

Frequency  response  constraints  allow  the  user  to  improve  such  system  characteristics 
as  command  following  and  stability.  The  frequency  constraints  are  applied  as  an 
upper  bound  •7ij(w)  to  the  magnitude  of  a  single  input  single  output  (SISO)  transfer 
function  Hij  in  the  multiple  input  multiple  output  (MIMO)  system  H: 

Ilf.iHI  <  7.iH  (2-28) 

where  Hij  is  parameterized  and  discretized  as  follows: 

Uu  JV— 1 

Hij{z)  —  Tijj  +  X/  53  '^2,ipFn<i>pqnTz,qj-  (2.29) 

p=l  qszl  n=:0 

Breaking  Hij{(jj)  down  into  real  and  imaginary  components,  the  constraint  in  Equa¬ 
tion  2.28  is  equivalent  to: 

Re[ifij(a;)]  cos  6  +  Im[^.v(w)]  sin 6  <  ^ij{io)  Me  €  [0, 27r).  (2.30) 

In  [5,  7],  this  form  of  magnitude  constraint  is  approximated  by  a  finite  number  of 
linear  constraints.  This  is  done  by  choosing  a  discrete  set  of  evenly  spaced  9  between 
0  and  27r: 

Re[Hij{u})]  cos  $„,N  +  Im[ffij(u;)]  sin  e„^N  <  7»i(^)  cos  ^  (2.31) 

where  ^  for  n  =  1, . . . ,  AT. 

After  defining  Sijpq{u;)  to  be: 

=  22,<p(‘*^)[Fo(w)  •  •  •  Fjv-i(w)]r3,gj(u;) 
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(2.32) 


and  defining  Lij{u)  to  be: 

=  7«i  ^  cos  9n,N  -  sin  ^h.n  (2.33) 

the  approximate  constraint  of  Equation  2.31  corresponds  to  a  finite  number  of  linear 
scalar  constraints  on  (f): 

Tiu 

cos  6n,N  +  Ini[‘S'.ip5(w)]  sin  6„,N)<l>pq  <  Lij{u}).  (2.34) 

p=i  g=i 

With  <f)  being  the  vector  of  <^pq’s,  Equation  2.34  is  now  written  as  Afreg<f>  <  &/reg  and 
can  be  applied  to  the  constrained  optimization  problem. 

2.4.2  Time  Response  Constraints 

Constraints  on  time  responses  allow  the  user  to  try  to  improve  output  characteristics 
to  a  given  input  disturbance.  For  a  closed  loop  MIMO  system  H{z),  time  constraints, 
in  the  form  of  upper  bounds  7„p  and  lower  bounds  'yiowi  applied  to  the  transient 
response  of  the  SISO  transfer  function  Hij{z): 

liow,ij[k]  <  Zij[k]  <  7„p,,i[fc].  (2.35) 

The  output  Zij  from  a  given  disturbance  wj  is  foimd  by  the  following  convolution: 

Zij[k]  =  (hij  *  u;j)[A:]  (2.36) 

where  is  the  impulse  response  of  Hij{z).  The  transient  response  constraint  of 

Equation  2.35  can  be  depicted  as  a  constraint  on  <f): 

Hu 

-  wj  *  ti,o)[fc]  <  £  £(«>J  *  rnijpq)[k]<f>pq  <  {7up,ij  -  Wj  *  ti,ij)[k]  (2.37) 

p=l9=l 

where  rriijpq  is  defined  in  Equation  2.16.  With  <t>  being  the  vector  of  <f>pq^s,  Equa¬ 
tion  2.37  is  now  written  as  Atime<f>  ^  bume  and  can  be  applied  to  the  constrained 
optimization  problem. 
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2.5  Model  Reduction 


As  discussed  in  [1],  a  post  optimization  controller  K  will  generally  have  the  order  of 
the  plant  plus  the  entire  order  of  Q,  which  is  based  on  the  number  of  basis  functions 
used  to  estimate  Q  in  Section  2.3.  Often,  the  resulting  K  has  too  high  an  order  for 
implementation.  If  the  user  does  not  wish  to  implement  the  high  order  K ,  a  form  of 
model  reduction  is  necessary  for  obtaining  a  reduced  compensator  Kred  of  acceptable 
order. 

There  are  many  forms  of  model  reduction.  For  the  purposes  of  this  thesis,  model 
reduction  will  be  reduced  to  two  options:  Balance  and  Truncate  and  Fractional  Bal¬ 
anced  Reduction  (FBR).  The  Balance  and  Truncate  method  is  chosen  for  its  widely 
accepted  use.  FBR,  an  extension  of  the  Balance  and  Truncate  method,  is  chosen 
because  of  its  applicability  to  models  with  unstable  modes. 

2.5.1  Balance  and  Truncate 

Balance  and  Truncate  is  the  most  widely  used  method  for  model  reduction.  Based 
on  Moore’s  algorithm  [13],  Balance  and  Truncate  relies  on  a  Cholesky  decomposition 
to  factor  the  controllability  and  observability  grammians,  which  then  balance  the 
system. 

Balanced  Realization 

The  first  step  of  Balance  and  Truncate  is  to  arrange  the  plant  into  an  internally 
balanced  form.  A  given  nth  order  minimal  stable  model  G(s)  is  defined  as  follows: 

G{s)  =  C{sl  -  A)-^B  (2.38) 

where  (A,  B,  C)  is  the  state  space  representation  of  G{s)  and  I  is  the  nth  order  identity 
matrix.  The  first  step  of  balancing  the  plant  is  to  solve  for  the  observability  grammian 
Wo  and  the  controllability  grammian  Wc  in  the  following  Lyapunov  functions: 

AWo+WoA^  +  BB'^  =  0 
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(2.39) 


A'^Wc  +  WcA  +  C'^C  =  0 


(2.40) 


where  Wo  and  Wc  are  unique,  symmetric,  and  positive  definite. 

The  next  step  is  to  decompose  Wo  with  the  Cholesky  decomposition: 

Wo  =  BF'R.  (2.41) 

This  is  used  to  form  the  matrix  RWcR^t  which  is  subsequently  diagonalized: 

RWcR'^  =  (2.42) 

where  U'^U  =  I  and  S  is  defined  as  a  diagonal  matrix  of  Hankel  singular  values, 
ordered  greatest  to  least  (ai  >  a2  >  . . .  >  <t„): 

S  =  diag{cri,cr2,...,<7„}.  (2.43) 

A  Balancing  Transformation  Matrix  ^bai  is  defined  as: 

$6a/  =  (2.44) 

Now  G(5)  can  be  transformed  into  the  minimal,  internally  balanced  model  Ghai{s)- 

Gbai(s)  =  Cm{sI  -  At^r'^Bbai  (2.45) 

where  Abai  =  ^baiA^^i,  Bbai  =  ^baiB,  and  Cbai  =  The  condition  for  any 

internally  balanced  plant,  such  as  Gbai{s),  is  that  the  corresponding  S  is  the  solution 
to  both  the  controllability  and  observabihty  Lyapunov  functions: 

A6a,S  +  SAL  +  BbaiBlt  =  AI{L  +  1:Am  +  ChCiai  =  0  (2.46) 

Truncation  of  Balanced  Model 

The  next  step  is  to  determine  the  size  of  the  reduced  model  r,  where  r  <n.  Now  the 
internally  balanced  Ghai{s)  can  be  represented  by: 

Abal  =  Bhal  =  Chal  =  [Gy  C'2]  (2-47) 

A21  A22  B2 
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where  An  is  size  r  x  r,  Bi  is  size  r  x  u  (u  inputs),  and  Ci  is  size  y  x  r  {y  outputs). 
Likewise,  E  is  also  partitioned: 


[  0  S2  J 


where  Si  =  diag{a-i, (T2,  • . . , (Xr)  and  S2  =  diag{ar+i, o-r+2)  •  •  •  > <x„}. 

Taking  the  components  of  Equation  2.47,  the  truncated  reduced  rth  order  model 

Gred  is  defined: 

Gred{s)  =  Ci{sl  -  An)-^Bi  (2.49) 

where  /  is  the  rth  order  identity  matrix.  This  tnmcated  rth  order  model  satisfies 
Equation  2.46  in  the  following  manner: 

AiiSi  +  El  A?;  +  BiBj  =  A[iSi  +  EjAn  +  C^Ci  =  0  (2.50) 

and  is,  therefore,  internally  balanced. 

Balance  and  TVuncate  only  reduces  stable  modes  of  the  model.  Unstable  modes 
are  removed  before  reduction  and  returned  to  the  model  after  reduction.  Balance 
and  Truncate  can  also  take  advantage  of  frequency  weighting  [5].  This  allows  the 
modes  to  be  emphasized  based  on  frequency.  If  a  lowpass  filter  is  used  as  a  frequency 
weight,  the  model  reduction  emphasizes,  and  likely  keeps,  the  low  frequency  states 
and  truncates  the  high  frequency  states.  Further  commentary  on  the  Balance  and 
Truncate  method  can  be  found  in  [5,  11,  12,  13]. 

2.5.2  Fractional  Balanced  Realization 

The  Fractional  Balanced  Reduction  method  of  model  reduction  was  developed  by 
Meyer  [11,  12],  extending  the  Balance  and  Truncate  method.  FBR  allows  for  reduc¬ 
tion  of  all  modes  of  a  system,  where  Balance  and  Truncate  fails  to  reduce  unstable 
modes.  FBR  takes  advantage  of  coprime  fractional  representations  of  the  plant,  gen¬ 
eralizing  Balance  and  Truncate  for  all  models. 

FBR  starts  with  the  minimal  state  space  representation  of  model  G(s),  as  de¬ 
scribed  in  Equation  2.38.  The  first  departure  from  Balance  and  Truncate  is  the  need 
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to  solve  the  following  algebraic  Riccati  equation  (ARE): 


PA  +  A'^P  -  PBB'^P  +  C^C  =  0  (2.51) 

where  P  is  the  positive  definite  solution.  This  step  is  based  on  a  coprime  fractional- 
ization,  with  the  derivation  found  in  [11,  12].  A  parameter  Ck  is  defined  as  follows: 

Ck  =  -B'^P.  (2.52) 

A  new  dynamics  matrix  A  is  defined  as: 

A  =  A  +  BCk.  (2.53) 


An  entirely  new  state  space  representation  is  then  defined  as: 

G{s)  =  C(sl  -  A)-^B 

1 

C 


(2.54) 


where  C  = 


Ck 


.  It  should  be  noted  that  Equation  2.51  is  now: 
PA  +  A^P  +  CICk  +  =  0 


(2.55) 


making  P  the  observability  grammian  of  this  system. 

FBR  now  returns  to  Balance  and  Truncate  procedures.  Using  the  procedure  in 
Section  2.5.1,  balance  the  system  G{s)  to  get  the  balanced  model  Gbai{s).  In  the 

Cbal 


balanced  model,  Chai  = 


C'ic, 


.6al 


Then  partition  Gbai(^)  ^.s  in  Equation  2.47, 


noting  that  Ci  = 


Cl 

Ck,\ 


T'he  second  departure  from  Balance  and  Truncate  is  the  partition  manipulation 
procedure,  based  on  coprime  fractions  and  Graph  Hankel  singular  vales,  derived  in  [11, 
12].  The  manipulation  is  simplified  to  defining  An  as  follows: 


All  =  .4.11  •"  BiCk,i 


(2.56) 


where  Bi  =  Bi.  Similar  to  Balance  and  Truncate,  the  rth  order  FBR  reduced  model 
Gred  is  now: 

Gred{s)  =  Ci{sI-An)-^Bi  (2.57) 

where  I  is  the  rth  order  identity  matrix. 
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Chapter  3 


Decoupled  Hierarchical  System 
Architecture 


This  chapter  presents  the  decoupled  hierarchical  system  architecture,  the  foundation 
of  the  solution  procedure  found  in  the  next  chapter.  Consider  the  example  system 
depicted  in  Figure  3.1,  there  are  four  channels,  one  for  each  Cartesian  position  (r,^,^) 
and  the  heading  Each  channel  has  an  associated  primary  control  actuator.  For 
example,  the  control  actuator  u\  is  associated  with  the  x  channel.  Since  each  channel 
has  an  associated  control  actuator,  it  is  assumed  each  channel  can  have  its  own 
associated  controller.  For  each  channel,  the  dynamics  are  entirely  coupled.  The 
first  decoupling  separates  the  dynamics  of  the  entire  system  into  individual  channels. 
Within  each  channel,  more  coupled  dynamics  exist.  For  the  purposes  of  this  thesis, 
a  coupled  system  will  be  the  coupled  channel  (which  is  decoupled  from  the  other 
channels)  as  depicted  in  Figure  3.2.  A  decoupled  hierarchical  system  will  be  decoupled 
within  the  channel.  Here,  the  pertinent  states  are  decoupled  from  one  another  as  in 
Figure  3.3. 

There  are  several  reasons  for  decoupling  the  system  into  the  hierarchical  structure. 
There  might  be  bandwidth  separation  between  the  modes  of  the  dynamics.  Decou¬ 
pling  would  give  the  control  engineer  the  advantage  of  control  design,  one  mode  at  a 
time,  with  the  opportunity  to  design  for  each  bandwidth  separately.  This  also  means 
that  each  mode  can  be  controlled  individually  using  multiple  control  techniques.  The 
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Figure  3.1:  Entire  System  of  Coupled  Channels 

hierarchical  structure  also  simplifies  the  problem.  Solving  several  decoupled  SISO 
(single  input/single  output)  systems  is  often  a  simpler  and  more  computationally 
efficient  process  than  solving  a  coupled  MIMO  system. 

3.1  System  Architecture 

A  coupled  MIMO  system  should  have  a  structure  similar  to  Figure  3.2.  In  this 
structure,  the  input  is  the  commanded  state  The  commanded  state  is  fed  into 
the  MIMO  plant  P{s),  where  the  coupled  states  (2:1,0:25^3)  are  the  outputs.  The 
outputs  are  fed  back  into  the  controller  K{s),  which  then  feeds  the  control  input  u 
into  the  plant.  Coupling  may  exist  between  internal  plant  states. 

For  the  purposes  of  this  thesis,  a  decoupled  system  will  have  the  hierarchical 
structure  depicted  in  Figure  3.3.  Here,  the  MIMO  system  is  broken  down  into  SISO 
inner  loops.  Each  pertinent  coupled  state  is  allocated  its  own  loop.  The  order  of  the 
loops  is  user  defined.  However,  it  is  suggested  that  the  state  loops  are  ordered  by 
having  the  faster  and  higher  order  derivative  states  in  the  inner  loops.  For  example, 
when  dealing  with  aircraft,  tracking  position  and  velocity  profiles  may  be  primary 
performance  objectives.  In  this  case,  the  faster  attitude  states  (roll,  pitch,  yaw, 
and  associated  rates)  are  allocated  to  the  more  internal  loops  and  the  position  states 
(position,  velocity,  and  acceleration)  are  allocated  to  the  outermost  loops.  In  addition, 
attitude  rate  state  loops  are  inside  the  attitude  state  loops  and  the  acceleration  loop 
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Figure  3.2:  Coupled  System  Architecture  for  a  Decoupled  Channel 

comes  before  the  velocity  loop,  which  comes  before  the  position  loop. 

The  plant  model  is  partitioned  into  a  series  of  transfer  functions  from  the  control 
u  to  the  pertinent  states.  The  error  of  the  innermost  state  sserr  is  the  input  for  the 
innermost  loop  controller  Ks^s),  that  determines  the  control  input  u.  The  control 
input  is  then  sent  to  the  innermost  plant  ^3(5),  which  sends  out  the  innermost  state 
as  an  input  to  the  next  loop’s  plant  P2{s)  and  as  negative  feedback  in  the  innermost 
loop.  In  the  subsequent  outer  loops,  the  commanded  states  become  the  outputs  of 
the  outer  controllers,  with  the  commanded  state  of  the  outer  loop  Xic  acting  as  the 
original  input  to  the  system.  Likewise,  the  outermost  state  xi  is  the  last  output  of 
the  system. 


3.2  Example  of  Coupled  Dynamics 

The  mathematical  decoupling  of  the  system  starts  with  the  state  space  description  of 
the  coupled  dynamics  P{s).  For  the  purposes  of  the  rest  of  this  chapter,  the  example 
system  will  be  the  system  depicted  in  Figure  3.2  and  Figure  3.3  with  three  pertinent 
states  {xi,X2,X3)  and  two  internal  states  (a:,-  and  xu).  The  state  space  representation 
of  the  coupled  system  P{s)  is: 

X  =  Ax  +  Bu 
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(3.1) 


Figure  3.3:  Decoupled  System  Architecture 


where  the  dynamics  matrix  A  and  control  matrix  B  are: 
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(3.2) 


with  the  control  as  a  sole  control  input  u  and  the  states  are  s  =  [xi  X3  Xi 


3.3  The  Innermost  Loop 

To  construct  the  plants  of  the  inner  loops,  the  dynamics  that  couple  the  pertinent 
states  are  removed  and  are  then  returned  with  each  outer  plant.  For  the  innermost 
loop,  the  system  is  truncated  into  the  innermost  pertinent  state  and  the  internal 

states.  The  input  to  the  system  is  the  commanded  innermost  state. 

■ 

The  generic  innermost  loop  is  depicted  in  Figure  3.4  for  R  pertinent  states.  =  3 
for  the  current  example,  where  the  commanded  state  input  is  The  state  space 
representation  of  the  innermost  system  ^3(5),  as  described  in  Figure  3.5,  has  the 
tnmcated  states  ^  =  [xs  Xi  This  state  space  system  is  as  follows: 


Xz  —  -^3^  +  .B3M3 

(3.3) 

+  D3.U3 

(3.4) 
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Figure  3.4:  Generic  Innermost  Loop 
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(3.6) 


with  the  inputs  —  [u  being  the  control  input  u  and.  the  commanded  state 

xzc-  The  outputs  =  [xserr  uY  axe  the  state  error  xserr  (which  is  fed  into  the 
controller),  the  pertinent  state  0:3,  and  the  control  u.  These  are  necessary  outputs  for 
applying  the  solution  procedure  and  analyzing  the  controller  design  using  the  SCTB 
in  Chapter  5. 


3.3.1  Design  of  Innermost  Controller  and  Closing  the  Loop 

Once  the  system  ^3(5)  is  defined,  the  innermost  controller  is  ready  for  design.  Fol¬ 
low  the  design  procedure  in  Chapter  4  to  solve  for  Kz-  The  state  space  representation 
of  Kz  is  as  follows: 

hzZ  =  -1-  BkzXZerr  =  -BkzXz  +  BkZ^Zc  +  AkS^z  (3-6) 

u  =  Ck3^z  +  DkzXZerr  =  —DkS^Z  +  DkS^Zc  +  Ck3^Z  (^•^) 

where  Xf.3  axe  the  controller  states  of  Kz  and  Akz  is  a,  rikz  X  rikz  matrix. 

After  Kz  is  designed,  close  the  loop  to  obtain  the  closed  loop  system  Hz'. 

Xjii  =  Auixiii  +  Biiixzc  (3*8) 
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Figure  3.5:  Innermost  Loop:  Model  Based  Representation 


=  CniXiii  +  DiijXac  (3.9) 

where 
Ain== 

(3.10) 

where  the  states  are  Xm  =  and  the  outputs  are  'gjjj  =  [xs  u]^.  The 

notation  >l3(:,  1)  corresponds  to  the  entire  first  column  of  matrix  As  and  A2{:,2  :  3) 
corresponds  to  a  matrix  consisting  of  the  second  through  third  coliunn  of  matrix 
^Ig.  The  input  is  the  commanded  state  xgc.  Once  the  closed  loop  formulation  of  the 
innermost  loop  is  completed,  the  next  loop  Gr-i{s)  (or  6^2(5)  for  the  example)  can 
be  formulated. 


As{:,l)  -  Bs{:,l)Dk3  A3{:,2:3)  Bg(:,l)C'fc3 

—Bks  0ix2  Aks 


Bill  = 


■®3(:>  l)L^ik3 


Cm  = 


1  0  0  0 
■Dks  0  0  Ck3 


Dni  = 


3.4  The  Intermediate  Loops 

For  each  intermediate  loop,  one  pertinent  state  is  added.  This  requires  that  state’s 
dynamics  to  be  returned  in  the  formulation  of  the  respective  intermediate  plant.  For 
the  intermediate  loop,  the  system  consists  of  the  previous  inner  pertinent  states  and 
the  internal  states.  The  input  to  the  system  is  the  commanded  intermediate  state. 
The  generic  intermediate  loop  is  depicted  in  Figure  3.6  for  pertinent  state  number 
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Figure  3.6:  Generic  Intermediate  Loop 


p.  The  current  example  has  p  =  2  and  X2c  as  the  commanded  intermediate  state.  The 
state  space  representation  of  the  innermost  system  ^2(5),  as  described  in  Figure  3.7, 
has  the  states  Xj  =  [x2  X3  Xi  xu  ^3]^.  This  state  space  system  is  as  follows: 


where 


A2  = 
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032 

Oi2 

0,i2 
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X2  =  A2X2  +  B2U2 
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(3.11) 

(3.12) 
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(3.13) 

with  the  inputs  U2  =  [xsc  ®2c]^  being  the  previous  commanded  state  xsc  and  the 
current  commanded  state  X2c>  The  outputs  =  {x2err  ®2  are  the  state  error 
iC2err  (which  is  fed  into  the  controller),  the  current  pertinent  state  X2,  and  the  control 
u.  As  explained  earlier,  these  outputs  are  chosen  for  application  of  the  solution 
procedure  and  analysis  of  the  controller  design  using  the  SCTB  in  Chapter  5. 
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Figure  3.7:  Intermediate  Loop:  Model  Based  Representation 


3.4.1  Design  of  Intermediate  Controller  and  Closing  the  Loop 

Once  the  intermediate  system  ^2(5)  is  defined,  the  intermediate  controller  K2{s)  is 
ready  for  design.  Follow  the  design  procedure  in  Chapter  4  to  solve  for  K2-  The  state 
space  representation  of  K2  is  as  follows: 


Xf.2  —  AhoXf.^  +  Bk2X2err  =  —Bk2X2  +  Bk2X2c  +  Ak2^2  (3-14) 

®3c  =  Ck2^2  +  Dk2X2err  =  —Dk2X2  +  Dk2^2c  +  Ck2Sjc2  (3.15) 


where  ^2  i'ii®  controller  states  of  ii'2* 

After  K2{s)  is  designed,  close  the  loop  to  obtain  the  closed  loop  system  H2‘ 


where 


An  = 


ill  =  AiiXii  +  BiiX2c 
Vll  =  CllXii  +  DllX2c 
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Figure  3.8:  Generic  Outer  Loop 

where  the  states  are  x//  =  [xj  ^2]^  outputs  are  =  [x2  u]'^.  The  input  is 

the  commanded  state  X2c.  Once  the  closed  loop  formulation  of  the  last  intermediate 
loop  is  completed,  the  outer  loop  ^1(5)  can  be  formulated. 

3.5  The  Outer  Loop 

For  the  outer  loop,  the  last  pertinent  state  Xi  is  added.  The  remaining  dynamics 
are  returned  in  the  formulation  of  the  outer  plant.  The  input  to  the  system  is  the 
commanded  outer  state  Xic. 

The  generic  outer  loop  is  depicted  in  Figure  3.6.  The  notation  is  the  same  for  the 
example.  The  state  space  representation  of  the  innermost  system  Gi(s)  as  described 
in  Figure  3.9  has  the  states  Xi  =  [x^  ^3  ^2]^-  This  state  space  system  is  as  follows: 

ki  =  AiXi  +  BiUi  (3.19) 


£1  =  CiXi  +  DiUi 


(3.20) 


Figure  3.9:  Outer  Loop:  Model  Based  Representation 
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(3.21) 

with  the  inputs  Ui  =  [x2c  ^ic]^  being  the  previous  commanded  state  X2c  the 
outer  commanded  state  xic-  The  outputs  =  [a:ierr  are  the  state  error 

^lerr  (which  is  fed  into  the  controller),  the  outer  pertinent  state  xi,  and  the  control 
u.  Once  again,  these  outputs  will  be  applied  to  the  solution  procedure  and  used  to 
analyze  the  controller  design  in  the  SCTB,  detailed  in  Chapter  5, 
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3.5.1  Design  of  Outer  Controller 

Once  the  outer  system  ^1(5)  is  defined,  the  final  controller  .^^1(5)  is  ready  for  design. 
Follow  the  design  procedure  in  Chapter  4  to  solve  for  K2.  This  should  be  all  the 
system  manipulation  necessary  for  the  decoupled  design.  Other  system  manipulation 
may  be  necessary  pending  model  reduction  implementation.  This  will  be  described 
in  Chapter  4. 
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Chapter  4 


Controller  Design  Solution 
Procedure 


A  flowchart  of  the  controller  design  solution  procedure  is  depicted  in  Figure  4.1. 
The  first  step  is  to  decouple  the  system  into  SISO  loops,  cts  described  in  Chapter  3. 
Once  the  innermost  decoupled  SISO  loop  is  defined,  as  in  Section  3.3,  then  that 
loop  is  ready  for  the  controller  design  solution  procedure  to  design  the  innermost 
controller.  Likewise,  for  subsequent  loops,  once  the  loop  is  defined  (in  Section  3.4  for 
an  intermediate  loop  or  Section  3.5  for  the  outer  loop),  the  controller  for  that  loop  is 
ready  for  the  design  solution  procedure. 

With  the  dynamics  set  up,  the  optimization  problem  is  ready  to  be  defined.  The 
first  step  in  defining  the  optimization  problem  is  to  choose  an  objective  function. 
Choosing  an  objective  function  is  described  in  detail  in  Section  4.1.  Once  the  objective 
function  is  determined,  the  constraints  are  applied.  The  application  of  constraints 
is  described  in  Section  4.2.  With  an  objective  function  ^d  constraints,  the  problem 
is  ready  for  the  Basis  Function  Algorithm,  which  solves  the  optimization  problem 
numerous  times  (with  a  different  basis  each  time)  to  find  the  controller  (the  design 
solution)  with  the  lowest  order.  This  algorithm  is  detailed  in  Section  4.3. 

After  the  Basis  Function  Algorithm  is  completed  for  a  set  of  constraints,  the  it¬ 
erative  portion  of  the  solution  procedure  begins.  If  the  Basis  Function  Algorithm 
determines  that  the  problem  is  infeasible  for  the  given  constraints,  then  the  con- 
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straints  need  to  be  relcixed  before  returning  to  the  Basis  Function  Algorithm.  If  the 
Basis  Function  Algorithm  has  a  feasible  solution,  then  the  solution  needs  to  be  ana¬ 
lyzed  for  its  performance  characteristics.  The  analysis  phase  is  coupled  with  model 
reduction.  Model  reduction  is  necessary  if  the  solution  is  not  of  implementable  order. 
Therefore,  the  original  solution  and  all  reduced  order  solutions  are  analyzed.  If  this 
analysis  determines  that  the  solution  (and  reduced  order  solutions)  has  unsatisfactory 
performance,  constraints  must  be  added  to  correct  this  performance  issue.  If  the  anal¬ 
ysis  determines  that  the  solution  (or  reduced  order  solution)  is  acceptable  but  could 
be  improved,  then  either  constraints  can  be  added  or  the  ^-minimization  objective 
function  may  be  attempted  in  conjunction  with  modest  constraint  additions. 

If  the  solution  controller  passes  the  scrutiny  of  the  analysis  step,  then  the  im- 
plementability  of  this  solution  needs  to  be  tested.  This  could  be  accomplished  on 
the  actual  hardware  the  model  used  for  the  design  is  based  on,  but  failures  on  the 
hardware  can  prove  costly.  Instead,  the  implementability  of  the  solution  should  be 
determined  on  a  simulator,  based  on  the  hardware.  If  the  simulator  test  proves  the 
solution  to  be  a  failure,  it  must  be  noted  why  and  how  the  solution  failed.  Then 
that  information  is  used  to  alter  the  constraints  for  another  iteration  of  the  solu¬ 
tion  procedure.  However,  if  the  solution  passes  the  simulator  test,  it  becomes  the 
controller  for  that  SISO  loop.  Now,  the  SISO  loop  is  closed  around  this  controller 
as  described  in  Chapter  3.  If  this  is  not  the  outermost  loop,  the  dynamics  of  the 
next  SISO  loop  are  added  and  prepared  for  the  impending  optimization  problem.  If 
the  loop  is  closed  aroimd  the  solution  controller  and  this  is  the  outermost  loop,  the 
procedure  is  completed. 

4.1  Choosing  an  Objective  Function 

Two  options  for  the  objective  function  are  the  'H2  objective  function  and  the  Q  min¬ 
imization  objective  function.  If  starting  from  scratch,  use  the  7^2  objective  function. 
If  modest  design  improvements  are  imposed  on  an  existing  controller  design,  the 
^-minimization  objective  function  may  be  appropriate.  However,  aggressive  design 
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Figure  4.2:  SISO  Loop  with  Low  Pass  Filter 
improvements  on  previous  designs  require  the  7^2  objective  function. 

4.1.1  'H2  Objective  Function 

Once  the  objective  function  is  chosen,  a  frequency  weight  should  be  applied  to  the 
SISO  loop.  This  solution  procedure  will  consistently  apply  the  frequency  weight  as 
a  first  order  lowpass  filter  (LPF)  on  the  commanded  input  Xc  as  in  Figure  4.2.  This 
is  to  counter  a  design  flaw  in  the  SCTB.  The  goal  of  the  optimization  is  to  minimize 
the  'H2  norm  of  the  transfer  function  from  the  commanded  state  to  the  error: 

The  SCTB  searches  for  a  characteristic  time  scale  within  its  solution  procedure.  In 
this  case,  there  is  unity  high  frequency  gain.  This  will  maJse  the  SCTB  act  as  if  there 
is  infinite  bandwidth  unless  there  is  a  LPF  on  the  input.  It  must  be  noted  that  the 
lowpass  filter  must  have  a  cutoff  frequency  of  less  than  half  of  the  sample  rate  of 
the  discrete  time  system.  For  example,  if  the  system  has  a  sample  rate  of  25  Hz, 
the  lowpass  filter  cutoff  frequency  must  be  less  than  12.5  Hz.  Implementation  of  the 
‘K2  objective  is  described  in  Section  2.2.1. 

4.1.2  Q-Minimization  Objective  Function 

The  ^-minimization  objective  may  be  used  only  if  there  is  a  very  good  nominal 
controller  design.  The  plant,  in  this  case,  must  be  the  closed  loop  system  Hnom  of 
the  original  plant  P  and  the  nominal  controller  Knom  S'S  in  Figure  2.4.  In  conjunction 
with  this  performance  objective,  constraints  need  to  reflect  slight  improvements  on  the 
nominal  design.  Constraints  that  reflect  more  than  modest  improvements  will  likely 
make  the  optimization  infeasible.  Implementation  of  the  Q-minimization  objective  is 
in  Section  2.2.2. 
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4.2  Setting  Constraints 


There  are  four  situations  that  lead  the  solution  procedure  to  the  constraint  setting 
phase.  The  initial  situation  immediately  follows  the  choice  of  objective  function. 
Here,  constraints  axe  applied  as  the  user  sees  fit,  usually  based  on  requirements.  The 
second  situation  will  follow  the  Basis  Function  Algorithm  if  it  determines  that  the 
optimization  problem  is  infeasible.  In  this  case,  it  is  assumed  the  problem  is  infeasible 
as  a  result  of  overly  aggressive  constraints.  Naturally,  constraints  must  be  relaxed 
for  the  next  iteration  to  try  to  get  a  feasible  solution.  The  third  situation  follows 
solution  controller  analysis.  If  the  analysis  determines  the  solution  is  unsatisfactory 
or  has  room  for  improvement,  constraints  need  to  be  altered  accordingly  for  the  next 
iteration.  The  final  situation  where  the  constraints  need  to  be  set  follows  the  simulator 
test.  If  the  solution  is  determined  to  be  either  unimplementable  or  lacking  robustness 
to  unmodeled  dynamics  and  sensor  noise,  constraints  must  be  altered  for  the  next 
iteration.  In  the  implementation  phase,  controllers  that  fail  usually  reveal  why  they 
failed,  whether  it  is  a  stability  failure  or  a  command  following  failure. 

The  types  of  constraints  available  for  this  solution  procedure  are  in  the  frequency 
domain  and  the  time  domain.  Frequency  domain  constraints  are  upper  bounds  on 
the  Bode  magnitude  for  a  given  input/output  pair.  Time  constraints  are  upper  and 
lower  bounds  on  the  response  curves  of  outputs  in  relation  to  a  given  input,  such  as 
a  unit  step  response. 

4.2.1  Frequency  Domain  Constraints 

Frequency  constraints  can  be  set  directly  to  the  desired  transfer  function.  The  de¬ 
sired  transfer  function  could  be  the  closed  loop  transfer  function  or  the  sensitivity 
transfer  function,  which  will  be  defined  in  terms  of  the  generic  SISO  loop  structure 
of  Figure  4.3.  For  the  purposes  of  this  procedure,  all  constraints  are  upper  bounds 
for  the  Bode  magnitude  of  the  given  transfer  function.  This  is  consistent  with  the 
SCTB.  Implementation  of  the  frequency  constraints  to  the  optimization  problem  is 
in  Section  2.4.1. 
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Figure  4.3:  Generic  SISO  Loop 
The  closed  loop  transfer  function  C(s)  is  defined  as: 

Constraining  the  high  frequency  end  of  the  closed  loop  transfer  function  should  im¬ 
prove  the  stability  robustness  of  the  system.  This  can  be  accomplished  with  con¬ 
straints  of  low  Bode  magnitude  on  the  high  frequency  region.  A  cluster  of  these 
constraints  can  create  a  roll  off  constraint. 

The  sensitivity  transfer  function  S(s)  is  defined  as: 

SM  =  ^  (4.2) 

Constraining  the  low  frequency  end  of  the  sensitivity  transfer  function  should  improve 
the  command  following  characteristics  of  the  system.  This  can  be  accomplished  with 
low  Bode  magnitude  constraints  on  the  low  frequency  region. 

4.2.2  Time  Domain  Constraints 

The,  time  constraints  are  applied  to  the  responses  of  the  desired  output  to  a  given 
input  disturbance.  Commonly,  the  given  input  disturbance  (or  commanded  state 
input)  will  be  in  the  form  of  an  impulse,  step,  or  ramp  response.  However,  it  is 
possible  to  have  a  user  defined  input. 

The  constraints  are  applied  as  upper  or  lower  bounds  to  the  specified  outputs, 
specifically  in  relation  to  the  chosen  input.  Typical  time  domain  constrains  are  on 
elements  such  as  maximum  response,  rise  time,  settling  time,  overshoot,  and  steady 
state  value.  Implementation  of  time  constraints  is  found  in  Section  2.4.2. 

In  this  thesis,  the  time  domain  responses  to  be  constrained  will  often  be  the  state 
error  curves,  with  a  desired  steady  state  value  of  zero.  The  input  response  of  choice 
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will  be  the  step  disturbance  response.  In  this  case,  the  commanded  state  input  is  set 
at  one  unit  value. 


4.3  Basis  Function  Algorithm  with  Optimization 
Solution 

The  Basis  Function  Algorithm  is  applied  for  each  set  of  desired  constraints.  This 
is  clearly  the  most  tedious  portion  of  the  solution  procedure.  The  main  goal  of  the 
Basis  Function  Algorithm  is  to  provide  a  feasible  solution  to  the  optimization  that 
results  in  the  lowest  order  controller.  To  achieve  this  goal  it  is  necessary  to  achieve 
some  subgoals. 

The  first  of  these  subgoals  is  to  reduce  the  number  of  basis  functions  used  to 
approximate  the  optimal  Q.  This  requires  a  good  basis  choice.  A  poor  basis  choice  can 
malce  the  problem  highly  complex  computationally,  requiring  an  inefficient  amount 
of  time  for  a  design.  In  addition,  the  problem  may  result  in  a  design  of  such  a  high 
order  that  the  controller  cannot  be  implemented,  rendering  the  optimization  useless. 
On  the  other  hand,  an  efficient  basis  will  result  in  lower  order  controllers  that  were 
designed  in  less  computational  time. 

Another  subgoal  is  based  on  the  notion  of  reducing  the  number  of  basis  functions, 
to  find  the  most  efficient  time  scale  for  those  forms  of  basis  that  depend  on  time 
scales.  As  described  in  Section  2.3,  there  are  many  forms  of  basis  functions.  The 
Legendre  and  Laguerre  functions  require  a  given  time  scale  to  define  the  basis.  The 
better  the  time  scale,  the  fewer  functions  necessary  to  approximate  the  optimal  Q. 
A  time  scale  is  better  than  another  if  it  requires  fewer  functions  to  achieve  a  feasible 
solution  or  it  uses  the  same  number  of  functions  with  a  better  objective  value. 

The  last  subgoal  is  to  reduce  the  amount  of  computational  time  necessary  for 
achieving  the  optimal  solution.  Finding  the  most  efficient  basis  requires  solving  the 
optimization  over  each  iteration  of  the  problem.  An  algorithm  is  necessary  to  make 
the  procedure  consistent  and  quicker  than  pure  brute  force. 
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4.3.1  Step  1:  Best  Homogeneous  Basis 

The  first  step  of  the  basis  function  algorithm  is  to  find  the  best  basis  using  only  one 
type  of  basis,  a  homogeneous  basis.  It  is  in  this  phase  of  the  algorithm  where  the  best 
time  steps  for  the  Legendre  and  Laguerre  basis  will  be  determined  for  future  steps. 

According  to  [5],  the  the  Legendre  basis  provides  the  best  initial  guess.  Using  that 
logic,  the  Basis  Function  Algorithm  will  use  the  Legendre  as  the  first  homogeneous 
basis.  Based  on  what  is  known  about  the  dynamics  of  the  loop  to  be  controlled,  the 
designer  must  make  an  educated  guess  at  a  time  scale  frequency  ato.  In  addition,  the 
user  must  have  an  idea  how  large  of  a  solution  will  be  tolerated.  An  upper  bound 
on  the  number  of  basis  functions  iVtd  is  implemented  for  this  reason.  Once  the  time 
scale  and  upper  bound  are  determined,  the  Algorithm  is  ready  to  enter  the  Basis 
Reduction  Procedure  (BRP). 

The  Basis  Reduction  Procedure 

The  BRP  applies  to  all  types  of  basis  functions,  not  just  the  Legendre.  The  homoge¬ 
neous  BRP  is  detailed  in  a  flowchart  in  Figure  4.4.  After  entering  the  BRP,  values  for 
the  current  number  of  basis  functions  N,  the  lower  limit  (the  Icirgest  known  infeasible 
number  of  basis  functions)  A,„/,  and  the  smallest  known  number  of  basis  functions 
that  provide  an  optimal  solution  (an  upper  limit)  Nbett  are  all  initialized.  iV,„/  is 
initialized  to  zero,  N^egt  is  initialized  to  infinity,  and  N  is  user  defined.  Beyond  the 
initial  iteration,  there  may  be  a  best  previous  solution,  say  to  a  previous  time  scale, 
that  used  Nsest  functions.  For  every  subsequent  iteration  entering  the  BRP,  Nbeat 
will  be  initialized  to  be  NBe$t  instead  of  infinity.  Since  Nsest  is  assumed  to  be  within 
Ntoi,  Nbett  would  then  replace  Ntoi  as  the  upper  limit  on  the  number  of  functions. 

After  initialization,  the  optimization  problem  is  solved  for  N  functions.  If  the 
feasible  solution  is  foimd,  then  Nbegt  will  be  set  to  be  N,  replacing  Ntoi  as  the  upper 
limit  on  N.  Then  N  is  decreased,  but  within  the  lower  limit  Ni„f.  If  the  problem  is 
infeasible  the  lower  bound  AT,,,/  is  raised  to  be  N.  Then  N  is  increased,  but  within  the 
upper  limit  Ntd  (or  Nbest)-  The  relationship  between  Ntd,  iV6c«e,  and  iV,„/  is  shown 
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Enter:  Must  know  objective, 
constraints,  basis  type,  Ntoi, 
and  time  scale  (if  applicable) 


Solve  for  N  functions 


Increase  N  Increase  N 
30  JV  <  Nbest  SO  N  <  Ntol 


N=Ninf\  \N  =  Ninf\ 


Feasible  Solution? 


\N  =  Nbest] 


N-X<Ni„,l 


N+X>Ntat^ 


Nbca  =  oo? 


Problem 

infeasible 


N  +  X>Nt^^ 


Solution  has  Nbest  functions 
(for  Nbest  =  Nsest,  best 
solution  is  previous  solution) 


Figure  4.4:  Flowchart  of  Homogeneous  Basis  Reduction  Procedure 
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Iteration  0  Iteration  1  Iteration  2  Iteration  3 


Figure  4.5:  Basis  Reduction  Procedure  Diagram 

in  Figure  4.5. 

This  recursion  continues,  closing  in  the  upper  and  lower  bounds  imtil  the  solution 
is  within  a  value  A  of  both  the  upper  and  lower  bounds.  The  more  precise  the  basis, 
the  lower  the  value  of  A.  However,  a  higher  A  can  reduce  the  computational  time. 
If  no  feasible  solutions  axe  fotmd  and  (N  +  A)  breahs  the  toleration  limit  Nt<j,  then 
the  problem  is  deemed  infeasible  and  the  constraints  need  to  be  relaxed.  If  a  feasible 
solution  is  found,  Ntest  =  Nsest  and  the  iteration  for  the  best  time  scale  can  begin. 

The  Time  Scale  Iteration 

The  iteration  for  the  frequency  of  the  best  Legendre  time  scale  is  similar  to  the  BRP. 
After  the  initial  solution  is  foimd,  find  the  smallest  pole  in  Q,  pq.  The  upper  and 
lower  bounds  on  the  time  scale  frequency  at  are  then: 

Pq  ^  ^Pq  (4-3) 

where  the  upper  bound  is  just  twice  the  lower  bound.  Each  new  value  of  at  that  results 
in  a  successful  solution  will  then  create  a  new  upper  or  lower  bound.  A  successful 
solution  can  be  determined  in  one  of  two  ways.  One  is  where  a  feasible  solution  is 
found  with  less  than  Nseat  functions  (and  then  Nteat  =  ^Beat)-  The  other  is  where  a 


60 


successful  solution  is  also  found  with  Nsett  functions,  but  the  objective  value  for  the 
optimization  is  lower  than  the  previous  Nsest  solution.  Like  the  BRP,  this  recursion 
continues,  closing  in  the  upper  and  lower  bounds  until  the  solution  is  within  a  value  Xa 
of  both  the  upper  and  lower  bounds.  The  time  scale  frequency  will  then  be  coronated 
as  aieg,  to  be  used  for  Legendre  functions  for  the  remainder  of  the  Basis  Function 
Algorithm.  The  more  precise  the  basis,  the  lower  the  value  of  Ao*  However,  a  higher 
Xa  can  reduce  the  computational  time. 

The  time  scale  iteration  for  the  Laguerre  basis  is  slightly  different  than  the  Leg¬ 
endre  basis.  There  are  no  set  upper  and  lower  bounds  on  the  time  scale  frequency,  as 
in  the  Legendre.  The  original  time  scale  frequency  ato  is  treated  as  the  upper  bound. 
Then  the  time  scale  frequency  is  systematically  decreased  tmtil  no  successful  solutions 
are  found.  The  lowest  time  scale  frequency  to  have  a  successful  solution  will  be  titled 
OLag-  If  decreasing  the  time  scale  frequency  below  ato  did  not  provide  a  successful 
solution,  then  ato  is  treated  as  the  lower  bound.  Then  the  time  scale  frequency  is 
systematically  increased  until  no  successful  solutions  are  found.  The  highest  time 
scale  frequency  to  have  a  successful  solution  will  be  titled  aiag-  If  increasing  the  time 
scale  frequency  above  ato  did  not  provide  a  successful  solution,  then  ato  is  coronated 
as  a  tag,  to  be  used  as  the  frequency  of  the  Laguerre  time  scale  for  the  remainder  of 
the  Basis  Function  Algorithm. 

Ordering  the  Basis  Types 

After  finding  the  best  homogeneous  Legendre  basis  solution,  the  best  homogeneous 
Laguerre  basis  is  determined.  The  Laguerre  basis  still  uses  the  BRP,  but  it  uses 
the  alterations  noted  in  the  time  scale  iteration.  The  order  of  the  solution  for  the 
Legendre  and  Laguerre  functions  are  compared  to  each  other.  Then  the  FIR  basis  is 
solved  for  one  more  than  the  amount  of  functions  used  in  the  higher  order  solution 
between  the  Legendre  and  Laguerre  basis.  For  example,  if  the  Legendre  basis  used 
11  functions  for  a  14th  order  solution  ^d  the  Legendre  basis  used  16  functions  for  a 
19th  order  solution,  then  the  FIR  basis  is  solved  for  17  functions  (one  more  than  the 
number  needed  for  the  higher  order  Laguerre  solution).  If  the  FIR  basis  produces  a 


feasible  solution,  which  it  rarely  does  in  this  case,  use  the  BRP  to  get  the  lowest  order 
FIR  solution.  K  the  FIR  basis  produces  an  infeasible  solution,  then  it  is  ordered  last 
out  of  the  three  types  of  functions. 

If  there  exists  a  good  original  controller  K„om  from  a  previous  iteration,  poles  from 
that  controller  may  be  added  or  form  the  basis  using  the  fixed  pole  model  (FPM) 
basis.  The  poles  are  ordered  from  smallest  to  greatest  (or  slowest  to  fastest).  The 
FPM  basis  is  a  special  case  and  is  not  involved  in  the  basis  ordering  process  unless  it 
is  the  most  efficient  basis. 

To  define  some  terms,  the  best  homogeneous  basis  form  will  be  called  Basis  Func¬ 
tion  A  creating  a  solution  controller  from  Na  functions  of  n^ith  order.  For  Legendre 
and  Laguerre  basis,  the  solution  controller  order  is  the  sum  of  the  plant  order  np,  the 
number  of  functions,  and  one: 

riA  =  Na  +  1  +  wp- 

For  the  FIR  basis,  the  solution  controller  order  is  the  sum  of  the  plant  order  and  the 
number  of  functions: 

nA  =  Na  +  np.  (4.5) 

The  next  best  homogeneous  basis  form  will  be  called  Basis  Function  B  and  the  worst 
homogeneous  basis  (most  likely  the  FIR)  will  be  called  Basis  Function  C. 

4.3.2  Step  2;  Best  Pair  Combination  Basis 

Once  the  homogeneous  solutions  axe  determined,  the  time  scales  are  finalized,  and 
the  basis  types  are  ordered  by  their  efficiency,  the  Basis  Function  Algorithm  is  ready 
to  find  the  best  combination  for  a  pair  of  basis  types.  The  goal  of  the  Pair  Combi¬ 
nation  Basis  Step  is  to  find  a  solution  lower  than  n>ith  order,  the  order  of  the  best 
homogeneous  solution.  A  successful  solution  will  find  a  feasible  solution  that  is  either 
lower  than  n>ith  order  or  is  n^th  order  with  a  lower  objective  function  value. 

Before  starting  this  step,  the  relationship  between  the  three  different  pair  combi¬ 
nations  must  be  understood.  For  Leg  Legendre  functions.  Lag  Laguerre  functions,  F 
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FIR  functions,  and  a  npth  order  plant,  the  controller  solution  order  n,  can  always  be 
determined  by  the  following  “2-combo’’  rules. 


n,  =  Lag  -F  Leg  -H  1  +  np 

(4.6) 

n,  =  Lag  -|-  F  -f  np 

(4.7) 

Ug  =  Leg  -|-  F  -|-  np 

(4.8) 

The  A/B  Combination 

The  step  starts  with  the  combination  of  basis  types  A  and  B.  Following  the  2-combo 
rules,  one  B-function  is  added  to  as  many  A-functions  necessary  Ao  to  comprise  the 
basis  that  keeps  n*  equal  to  If  this  basis  makes  the  optimization  infeasible,  then 
a  basis  of  two  B-functions  and  (Aq  —  1)  A-functions  is  tried.  If  the  optimization  is 
still  infeasible,  continue  this  pattern  with  a  basis  of  three  B-functions  and  (Aq  —  2) 
A-functions.  If  all  attempts  up  to  this  point  are  infeasible,  experience  has  indicated 
that  it  is  unlikely  to  find  a  feasible  solution  continuing  on  this  pattern.  Instead,  go 
straight  to  the  Alternate  Pair  Combinations  portion  of  this  section. 

If  the  basis  of  one  B-function  and  Ao  A-functions  provides  a  successful  solution, 
then  try  the  basis  of  one  B-function  and  (Ao  ~  1)  A-functions.  Keep  reducing  the 
amount  of  A-functions  by  one  until  the  optimization  is  infeasible.  The  best  solution 
from  this  iteration  is  n^ioth  order.  Following  the  2-combo  rules,  the  two  B-function 
basis  is  attempted  where  Ai  A-functions  are  used  such  that  n,  equals  riAo-  The 
same  procedure  applies  where  the  amount  of  A-functions  is  reduced  by  one  until  the 
optimization  is  infeasible.  This  general  procedure  is  applied  tmtil  the  addition  of  two 
B-functions  provides  no  successful  solutions.  At  this  point,  try  the  Alternate  Pair 
Combinations. 

Alternate  Pair  Combinations 

After  trying  the  conventional  combinations  for  pairs  of  basis  types,  these  two  alternate 
attempts  are  worth  trying.  First,  try  flipping  the  amotmt  of  the  superior  and  inferior 
functions.  In  the  A/B  case,  start  with  one  A-function  and  Aq  B-functions  (or  however 
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many  B-functions  necessary  to  match  the  order  of  the  current  best  solution).  Then 
follow  the  procedure  as  in  the  A/B  Combination  portion  of  this  section.  If  this 
provides  no  successful  solution,  only  two  iterations  were  wasted. 

The  other  alternate  attempt  is  to  try  an  even  number  of  A-functions  and  B- 
functions.  If  very  small  numbers  of  basis  functions  are  being  attempted,  this  might 
have  already  been  tried.  However,  for  large  numbers  of  basis  functions,  this  gives  the 
design  engineer  the  ability  to  check  an  A-heavy  basis,  a  B-heavy  basis,  and  now  a 
A/B  balanced  basis.  If  this  produces  a  successful  solution,  then  follow  the  procedure 
as  described  in  the  A/B  Combination  portion  of  this  section.  However,  if  this  does 
not  produce  a  successful  solution,  end  the  procedure  for  this  particular  basis  pair 
combination. 

A/C  and  B/C  Combinations 

The  A/C  and  B/C  combinations  are  foimd  in  the  same  manner  as  the  A/B  combi¬ 
nation,  with  the  C-function  assuming  the  inferior  function  role  (B  was  the  inferior 
function  in  the  A/B  Combination).  The  goal  for  these  combinations  is  finding  suc¬ 
cessful  solutions  that  improve  on  the  best  current  solution.  At  this  point,  the  best 
current  solution  is  either  the  n^ith  order  homogeneous  A-basis  solution  or  the  A/B  so¬ 
lution  that  improved  upon  the  n^th  order  solution.  The  Alternate  Pair  Combinations 
also  apply  to  these  combinations. 

Adding  Fixed  Pole  Model  Functions 

Using  FPM  functions  from  a  good  controller  from  a  previous  iteration  is  a  special 
Pair  Combinations  case.  Unless  the  FPM  functions  were  the  best  homogeneous  basis, 
these  functions  will  always  be  treated  as  the  inferior  function.  First,  the  A/FPM 
combination  is  tried  in  the  same  manner  as  the  A/B  combination.  If  no  successful 
solution  comes  out  of  an  A/FPM  combination,  then  the  Pair  Combination  procedure 
ends  there.  If  a  successful  A/FPM  solution  is  found,  then  the  B/FPM  combination 
is  tried.  The  same  conditions  apply  for  the  C/FPM  combination  to  be  attempted. 


4.3.3  Step  3:  Best  Three-Combination  Basis 

After  the  Pair  Combination  procedure  is  complete,  the  current  best  solution  will 
either  be  the  A-function  homogeneous  case  or  the  best  Pair  Combination  solution. 
The  best  Three-Combination  procedure  tries  to  determine  if  a  basis  combination  of 
Legendre,  Laguerre,  and  FIR  functions  can  improve  on  the  best  current  solution. 
If  the  Pair  Combination  is  the  best  current  solution,  the  procedure  is  simple.  Add 
one  Third  Function  and  remove  one  function  from  the  most  represented  of  the  Pair 
Combination.  If  this  produces  a  successful  solution,  a  function  of  the  most  represented 
Pair  Combination  is  removed  one  at  a  time.  This  continues  until  the  problem  is 
infeasible.  Then  another  Third  Function  is  added  to  the  basis  and  the  most  numerous 
Pair  Combination  member  loses  a  function.  This  procedure  continues  until  adding  2 
Third  Functions  does  not  produce  a  successful  solution. 

If  the  homogeneous  case  is  the  current  best  solution,  then  start  with  one  B- 
function,  one  C-function,  and  eis  many  A-functions  necessary  to  get  an  nj4th  order 
solution.  When  there  is  a  successful  solution,  remove  one  A-function  at  a  time  until 
the  problem  is  infeasible.  When  there  is  no  successful  solution,  remove  one  A-function 
and  add  one  B-function.  If  this  does  not  provide  a  successful  solution,  remove  the 
new  B-function  and  add  one  C-function  instead.  The  procedure  ends  when  two  B- 
functions  or  C-functions  are  added  without  producing  a  successful  solution. 

Adding  FPM 

For  a  known  given  good  controller,  if  no  Three-Combination  produces  a  solution,  try 
the  Three-Combination  with  the  FPM  functions  replacing  the  C-functions.  If  there 
is  a  Three- Combination  solution,  FPM  functions  can  be  added  as  a  fourth  form  of 
basis.  This  is  done  by  removing  the  one  function  of  the  most  represented  function 
from  the  Three-Combination  basis  for  every  FPM  function  added  to  the  basis.  When 
solutions  axe  found,  the  most  represented  function  reduces  one  at  a  time  until  the 
problem  is  infeasible.  When  adding  two  FPM  functions  fails  to  produce  a  successful 
solution,  the  procedure  is  over. 
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4.4  Analysis  and  Model  Reduction 


After  the  best  feasible  solution  is  found,  it  must  be  determined  whether  this  design 
truly  meets  all  the  constraints.  There  also  may  be  characteristics  of  the  design  that 
were  unaccounted  for  by  the  optimization.  In  addition,  there  might  be  an  imple¬ 
mentation  requirement  on  the  order  of  the  controller.  Very  often,  the  optimization 
will  produce,  at  best,  a  design  that  is  still  of  very  high  order.  This  is  where  model 
reduction  is  applied.  Model  reduction  is  also  useful  for  designs  of  moderate  order, 
because  it  is  always  better  to  have  a  lower  order  design. 

The  analysis  portion  of  the  procedure  determines  if  the  solution’s  performance 
is  acceptable.  If  the  solution  has  satisfactory  performance,  then  model  reduction 
is  applied  to  the  solution.  If  the  controller  has  all  stable  modes,  then  the  Balance 
and  Truncate  method  of  Section  2.5.1  will  reduce  the  system.  If  the  controller  has 
unstable  states  or  modes,  then  the  reduction  method  will  be  the  Fractional  Balanced 
Reduction  (FBR)  method  of  Section  2.5.2.  The  FBR  method  is  an  extension  of  the 
Balance  and  Truncate  method  for  unstable  systems.  Once  the  controller  is  reduced  to 
a  desired  order,  this  reduced  order  design  is  analyzed.  If  there  is  an  implementation 
requirement  having  controllers  less  than  n,th  order,  then  all  controllers  under  n,th 
order,  that  have  acceptable  performance,  advance  to  the  simulator  test  phase.  If  there 
are  no  designs  that  meet  this  requirement,  constraints  need  to  be  altered  for  a  new 
iteration  of  controller  design. 
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Chapter  5 


Structured  Control  Toolbox 


The  Charles  Stark  Draper  Laboratory  developed  a  MATLAB  based  Computer-Aided 
Design  (CAD)  tool  called  the  Structural  Control  Toolbox  (SCTB)  [5].  The  SCTB 
is  driven  in  a  Graphical  User  Interface  (GUI)  environment  that  greatly  reduces  the 
interaction  between  the  user  and  the  MATLAB  command  line.  Although  the  SCTB 
has  numerous  capabilities,  this  chapter  will  only  discuss  those  capabilities  necessary 
for  implementing  the  solution  procedure  presented  in  the  previous  chapter.  Additional 
details  and  alternate  explanations  of  the  SCTB  are  found  in  [5,  7]. 

5.1  Main  Panel 

Once  the  SCTB  is  launched  from  the  MATLAB  command  line,  an  empty  Main  Panel 
is  displayed,  as  in  Figure  5.1.  The  Main  Panel  has  several  menus.  On  the  top  menu 
bar,  there  are  four  pull-down  menus:  File,  Windows,  Options,  and  Help.  Within 
the  Main  Panel  there  are  pull-down  menu  bars:  Modify  (for  plants  and  controllers), 
Pre-Processor  (for  plants  and  controllers),  Design  Methods,  and  Analysis.  In 
addition  to  menu  bars,  there  are  function  buttons  on  the  Main  Panel.  The  function 
buttons  of  interest  include:  Import  P,  Close  Loop,  and  Close. 
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Figure  5.1:  SCTB  Main  Panel 

5.1.1  Top  Pull-Down  Menus 

The  File  pull-down  menu  has  six  items:  New,  Open  DB,  Close,  Print,  Save, 
and  Save  As.  The  New  option  will  clear  the  current  database  from  the  SCTB, 
giving  the  user  the  option  of  saving  or  deleting  the  current  database.  The  Open  DB 
option  will  load  a  previously  saved  database  into  the  SCTB.  If  there  is  already  an 
active  database,  the  Open  DB  command  will  replace  the  current  database,  giving 
the  user  the  option  of  saving  or  deleting  that  database.  The  Close  option  will  quit 
the  SCTB,  giving  the  user  the  option  of  saving  or  deleting  the  current  database.  The 
Print  option  will  print  the  panel  to  a  printer  or  file.  The  Save  and  Save  As  options 
will  save  the  current  database  to  a  MATLAB  MAT  file.  Under  the  Save  option,  the 
user  may  save  the  database  under  the  current  filename  or  save  the  database  to  a  new 
filename.  The  Save  As  option  saves  the  database  to  a  new  filename. 

The  Windows  pull-down  menu  lists  the  currently  opened  SCTB  figure  windows. 
Selecting  a  window  will  activate  that  particular  figure  window.  The  Options  pull¬ 
down  menu  has  two  items:  Empty  Trash  Can  and  Refresh.  The  Empty  Ttash 
Can  option  deletes  all  data  in  the  trash  can,  located  below  the  Analysis  menu  bar 
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in  the  Main  Panel.  The  Refresh  option  updates  the  Main  Panel  based  on  data  in 
the  current  database.  The  Help  pull-down  menu  lists  different  help  windows  that 
will  display  help  topics  on  the  SCTB,  MATLAB,  and  the  Main  Panel  itself. 

5.1.2  Panel  Menu  Bars 

The  Modify  menu  bars  have  three  items:  Group,  Rdit,  and  Modify  via  Simulink. 
The  Edit  option  is  the  only  option  of  interest  to  this  thesis.  The  Edit  option  allows 
the  user  to  modify  properties  of  the  selected  plant  or  controller  model  in  the  current 
database  using  portions  of  the  Data  Input  Module,  described  in  Section  5.2.  The 
Pre-Processor  menu  bars  have  two  items:  Reduction  and  Unity  Feedback.  The 
Reduction  option  is  the  only  option  of  interest  to  this  thesis.  The  Reduction  option 
allows  the  user  to  reduce  the  order  of  the  selected  plant  or  controller  model  in  the 
current  database  using  the  model  reduction  tool,  described  in  Section  5.5. 

The  Design  Methods  menu  bar  has  five  items:  Classical  Design,  Constrained 
Optim,  H2  Design,  Hinf  Design,  and  Mu  Tools.  The  only  option  of  interest  to 
this  thesis  is  Constrained  Optim.  This  option,  detailed  in  Section  5.3  allows  the 
user  to  design  a  controller  using  constrained  convex  optimization  techniques.  The 
other  options  are  all  various  other  forms  of  controller  design.  The  Analysis  menu 
bar  has  two  items:  Time/iVeq  Analysis  and  Signal  TYansmission  Analysis. 
The  Time/Freq  Analysis  item  is  the  only  option  of  interest  to  this  thesis  and  is 
detailed  in  Section  5.4.  This  allows  the  user  to  analyze  the  current  database  plants 
and  controllers  in  the  time  and  frequency  domains. 

5.1.3  Main  Panel  Buttons 

There  are  six  Main  Panel  buttons:  Import  P,  Import  K,  Close  Loop,  Real-Time 
K,  Solve  K,  and  Close.  The  buttons  of  interest  are  Import  P,  Close  Loop,  and 
Close.  The  Import  P  button  allows  the  user  to  load  a  plant  model  from  a  MATLAB 
MAT  file  into  the  SCTB  database  using  the  Data  Input  Module.  This  loading  process 
is  presented  in  Section  5.2.  The  Close  Loop  button  will  generate  the  closed  loop 
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Figure  5.2:  Data  Input  Module:  MAT  File  Selection 


plant  model  for  the  active  plant  and  the  active  controller  using  the  Closed  Loop  Tool, 
presented  in  Section  5.6.  The  Close  button  has  the  same  effect  as  selecting  the  Close 
item  of  the  File  pull-down  menu:  quitting  the  SCTB  and  giving  the  user  the  option 
of  saving  or  deleting  the  current  database. 


5.2  Loading  a  Plant  Model 

A  state  space  representation  of  a  plant  model  can  be  loaded  from  a  MATLAB  MAT 

■  .H-  • 

file  and  placed  into  the  SCTB.  When  the  Import  P  button  is  selected  in  the  Main 
Panel  the  first  Data  Input  Module  GUI  is  displayed  as  in  Figure  5.2.  Here,  the  GUI 
calls  for  the  MAT  file  of  the  plant  model.  Upon  receiving  this  file,  the  GUI  adds  the 
Select  the  type  of  data  portion  to  the  window  as  in  Figure  5.3.  This  portion  of  the 
GUI  has  three  buttons:  State-Space  model,  FRF  data,  and  ZPK  model.  For 
this  thesis,  only  state  space  models  are  used  and  the  State-Space  model  is  the  only 
button  of  interest  and  should  be  selected.  However,  if  the  model  was  in  frequency 
response  function  (FltF)  or  zero/pole/gain  (ZPK)  form,  the  other  two  options  should 
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Figure  5.3:  Data  Input  Module:  Data  Type  Selection 


be  selected. 

Once  the  State-Space  model  button  is  selected,  the  State-Space  Model  Input 
Window  is  opened  as  in  Figure  5.4.  The  variables  in  the  MAT  file  are  assigned  to  the 
variables  A,  B,  C,  D,  and  T,  in  Figure  5.4  to  build  the  plant  in  the  SCTB  database. 
The  A,  B,  C,  and  D  variables  correspond  to  state  space  matrices  and  T,  is  the  time 
step  for  discrete  systems.  For  continuous  models,  a  0  is  entered  in  the  Tg  field.  If  the 
state  space  matrices  in  the  MAT  file  are  not  named  A  ot  a,  B  or  6,  C  or  c,  and  D 
or  d,  respectfully,  the  SCTB  will  not  place  the  matrices  in  the  appropriate  field  and 
the  user  will  need  to  t3^e  the  proper  variable  names  into  their  corresponding  matrix 
fields.  Based  on  the  matrices  in  the  A,  B,  C,  and  D  fields,  the  corresponding  inputs 
and  outputs  are  listed  in  the  Input /Output  Descriptions  portion  of  the  GUI.  There 
the  user  may  edit  the  name  of  ^y  input  or  output.  Once  the  variables  are  set  and  the 
inputs  and  outputs  are  named,  the  Build  P  button  shordd  be  selected  to  continue 
assembling  the  plant  model. 


Figure  5.4:  Data  Input  Module:  State-Space  Model  Input  Window 

5.2.1  Defining  Plant  Characteristics 

The  Build  P  button  will  invoke  the  Plant  Definition  Window  as  in  Figure  5.5.  This 
window  allows  the  user  to  select  input  and  output  channels,  add  plant  uncertainty 
and  input /output  augmentation,  choose  between  discrete  and  continuous  time,  and 
edit  the  plant  name.  The  input  channel  selection  portion  of  this  panel  assigns  each 
input  to  be  a  control  actuator  (u),  a  disturbance  (d),  or  both.  At  least  one  input 
must  be  designated  a  control  actuator.  The  output  channel  selection  portion  of  this 
panel  assigns  each  output  to  be  a  sensor  (y),  a  measurement  (e),  or  both.  At  least  one 
output  must  be  designated  a  sensor.  This  thesis  will  not  select  any  plant  uncertainty, 
but  will  always  select  the  Add  Noise  Inputs  and  Add  Control  Outputs  buttons. 
Constrained  optimization  problems  involve  discrete  time  systems.  For  consistency, 
the  Tustin  conversion  method  is  always  selected  when  converting  continuous  systems 
to  discrete  systems.  Upon  completion,  selecting  the  Done  button  will  create  the 
plant  model,  add  it  to  the  SCTB  database  and  Main  Panel,  and  close  the  Data  Input 
Module  GUI. 
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Figure  5.5:  Data  Input  Module:  Plant  Definition  Window 

5.3  Constrained  Optimization  Controller  Design 

After  the  plant  is  defined  and  added  to  the  SCTB  database,  it  is  time  for  controller 
design.  Select  the  Constrained  Optim  item  of  the  Design  Methods  menu  bar  in 
the  main  panel. 

This  invokes  the  Constrained  Optimization  Controller  Design  panel  in  Figure  5.6. 
This  panel  has  four  main  modes:  Set  Objective,  Edit  Basis,  Edit  Freq.  Const., 
and  Edit  Time  Const.  There  axe  also  four  buttons  (excluding  the  Close  button): 
Feedback  Sign:+,  View  Controller,  SaVe,  and  Solve.  The  feedback  sign  toggle 
switch  should  always  remain  on  positive.  The  View  Controller  button  gives  you  a 
graphical  representation  of  a  solved  controller.  The  Save  button  saves  the  current 
solved  controller.  The  Solve  button  solves  the  optimization  for  the  given  objective, 
basis,  and  constraints. 
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Figure  5.6:  Constrained  Optimization  Controller  Design 

5.3.1  Setting  the  Objective  Function 

Upon  its  initial  opening,  the  Constrained  Optimization  Controller  Design  is  in  Set 
Objective  mode.  There  is  a  Norm  Optimization  menu  bar  requesting  the  user  to 
Select  a  Method.  The  menu  bar  has  three  items:  H2  Minimization,  Q  Minimiza¬ 
tion,  and  LI  Minimization.  The  7^2  ^-Minimization  Objective  Functions  are 
the  only  objectives  of  interest  to  this  thesis. 

7^2  .Minimization 

After  selecting  the  H2  Minimization  item  in  the  Norm  Optimization  menu  bar,  the 
window  looks  like  Figure  5.7.  Here,  the  inputs  and  outputs  are  based  on  the  generic 
SISO  loop  of  Figure  4.3.  There  are  two  parts  to  preparing  the  7^2  Objective  Function: 
selecting  the  input/output  pairs  and  adding  the  frequency  weights.  Selecting  the 
input/output  pairs  are  relatively  eeisy.  First,  the  user  should  select  the  input /output 
pair  that  needs  to  be  minimized.  In  the  generic  SISO  case,  x^rr  is  to  be  minimized, 
so  the  (xc,Xerr)  paif  shorJd  be  selected.  Also,  the  plant  is  assumed  to  be  four  block, 
meaning  that  the  number  of  regulated  outputs  exceeds  the  number  of  measurements 


Figure  5.7:  %2  Objective  Function  Control  Panel 

and  the  number  of  exogenous  inputs  exceeds  the  number  of  controls.  Therefore,  the 
pair  of  the  regulated  output  x„rNOISE  and  the  exogenous  input  uCONTROL  is 
also  selected.  The  intersections  of  the  two  selected  pairs  are  automatically  selected. 
This  is  reflected  in  Figure  5.7. 

Frequency  Weights 

To  add  frequency  weights,  the  FVequency  Weight  button  must  be  selected.  This 
will  invoke  the  Defining  Weighting  Matrices  GUI.  To  get  a  LPF  on  the  commanded 
input,  as  in  Figure  4.2,  select  the  W3  (Disturbance  Input-d)  button.  Then  the 
Fdit  button  will  invoke  the  Function  Editing  window.  Select  Lowpass  in  the  Filter 
menu  and  Order  (1  or  2)  in  the  Filter  Parameters  menu.  To  get  a  first  order 
LPF,  select  1st  Order  in  the  menu  bar.  To  add  the  desired  break  frequency,  choose 
Frequency  in  the  Filter  Parameters  menu.  Select  the  Update  button  to  enter  the 
LPF  into  the  General  Description  window.  Now,  select  the  Save  button  to  enter 
the  LPF  as  the  W3  function  and  return  to  the  Defining  Weighting  Ma,trices  panel. 
The  Bode  plot  of  the  LPF  should  be  displayed.  Selecting  the  Save  button  will  add 
the  frequency  weighting  to  the  objective  function  and  return  to  the  Objective 
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Figure  5.8:  Basis  Function  Control  Panel 
Function  Control  Panel. 

^-Minimization 

If  Q-minimization  is  the  desired  objective  function,  the  plant  must  be  a  product  of  the 
Close  Loop  Tool  in  the  Main  Panel.  The  Close  Loop  Tool,  presented  in  Section  5.6, 
will  produce  a  closed  loop  plant  consisting  of  an  open  loop  plant  and  a  nominal 
controller.  After  the  closed  loop  plant  is  selected  and  the  Constrained  Optimiza¬ 
tion-Controller  Design  Window  is  opened,  the  ^-Minimization  objective  is  chosen  by 
choosing  the  Q  Minimization  item  in  the  Norm  Optimization  menu  bar  of  the  Set 
Objective  panel. 

5.3.2  Editing  the  Basis 

When  the  objective  function  is  determined,  the  basis  may  be  edited  by  selecting  the 
Edit  Basis  button.  This  will  invoke  the  Basis  Function  Control  Panel  in  Figure  5.8. 
To  create  a  basis,  basis  functions  must  be  chosen.  If  Laguerre  or  Legendre  functions 
axe  desired,  a  time  scale  frequency  must  be  defined.  To  define  a  time  scale,  select  the 
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Figure  5.9:  Frequency  Constraint  Transfer  Function  Selection  Panel 

respective  Add  Frequency  (Hz)  button.  Tbe  appropriate  frequency  can  now  be 
entered.  Then  to  add  Laguerre,  Legendre,  or  FIR  functions,  enter  the  desired  number 
of  functions  in  the  respective  portions  of  the  panel.  To  enter  the  user  defined  FPM 
functions,  they  can  be  tahen  from  the  plant  (select  the  FVom  Plant  button),  from  the 
Q  of  the  last  solution  (select  the  R*om  Q  button),  or  from  the  workspace.  There  will 
be  a  list  of  pole  choices  in  the  Select  Poles  menu.  Select  the  desired  poles  and  add 
them  to  the  basis  by  selecting  the  Add  ==>  button.  Once  all  of  the  FIR,  Laguerre, 
Legendre  and  FPM  functions  axe  determined,  the  basis  is  created  by  selecting  the 
Create  Basis  State-Space  Model  button. 

5.3.3  Editing  Frequency  Constraints 

To  set  frequency  constraints,  select  the  Edit  IVeq.  Const,  button.  This  will  invoke 
the  Frequency  Constraint  Transfer  Function  Selection  panel  in  Figure  5.9.  To  select 
the  sensitivity  transfer  function  select  the  (xcXcrr)  pair.  To  select  the  closed 

loop  transfer  function  select  the  {xc,x)  pair.  Once  a  transfer  function  is  chosen, 
selecting  the  Add/Edit  Constraints  button  will  invoke  the  Frequency  Constraint 
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Figure  5,10:  Frequency  Constraint  Tool 

Tool  in  Figure  5.10.  The  Frequency  Constraint  Tool  allows  the  user  to  place  upper 
bound  constraints  on  the  Bode  magnitude  of  the  chosen  transfer  function.  These 
constraints  can  be  a  point  or  a  line  of  points.  The  point  density  of  the  constraint 
line  is  offered  in  terms  of  points  per  decade  and  is  chosen  in  the  constraint  portion  of 
the  tool,  just  under  the  Add  Points  and  Add  Lines  buttons.  After  the  constraint 
type  (point  or  line)  is  determined,  the  constraints  are  set  by  placing  them  on  the 
magnitude  grid  itself.  The  Tool  also  allows  the  user  to  plot  the  open  and  closed 
loop  versions  of  the  selected  transfer  function  to  compaxe  with  the  constraints.  In 
Figure  5.10,  a  high  frequency  roll  off  constraint  is  placed  on  the  closed  loop  transfer 
function 

5.3.4  Editing  Time  Constraints 

To  set  constraints  on  the  time  response,  select  the  £}dit  Time  Const,  button.  This 
will  invoke  the  Time  Constraint  Transfer  Function  Selection  panel  in  Figure  5.11.  To 
select  the  error  response  to  the  commanded  input  disturbance,  select  the  (ajcj^err)  pair. 
Once  a  transfer  function  is  chosen,  selecting  the  Add/Edit  Constraints  button  will 
invoke  the  Time  Constraint  Tool  in  Figure  6.12.  The  Time  Constraint  Tool  allows 
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Figure  5.11:  Time  Constraint  Transfer  Function  Selection  Panel 


Input  i^lOiUputJ^  Leaond 


Figure  5.12:  Time  Constraint  Tool 
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the  user  to  place  upper  and  lower  bounds  on  the  time  response  of  an  output  to 
a  given  input  disturbance.  The  Time  Constraint  Tool  has  three  menu  bars.  The 
Figure  menu  bar  toggles  between  the  constraint  plot  and  the  disturbance  plot.  The 
Disturbance  Type  menu  bar  allows  the  user  to  choose  the  type  of  input  disturbance. 
The  choices  are  Impulse,  Step,  Ramp,  and  Other.  The  Other  choice  is  a  user 
defined  input  disturbance.  The  third  menu  bar  toggles  between  Upper  and  Lower 
Constraints. 

These  constraints  can  be  a  point  or  a  line  of  points.  The  point  density  of  the 
constraint  line  is  in  terms  of  time  steps  per  constraint  (for  discrete  systems)  and 
is  chosen  in  the  Edit  Constraints  portion  of  the  tool.  After  a  disturbance  type  is 
chosen  and  a  constraint  type  (point  or  line)  is  determined,  the  constraints  are  set 
by  placing  them  on  the  magnitude  grid  itself.  The  Tool  also  allows  the  user  to  plot 
the  open  and  closed  loop  transients  to  compare  with  the  constraints.  In  Figure  5.12, 
time  constraints  are  placed  on  the  error  response  Xerr{t)  to  a  unit  step  commanded 
input  ®c(t)-  The  constraints  are  a  ±1  amplitude  constraint  and  a  settling  constraint 
of  ±0.01  from  8  seconds  to  15  seconds. 

5.3.5  Solving  the  Optimization  Problem 

After  the  objective,  basis,  and  constraints  have  been  set,  the  optimization  problem 
is  ready  to  be  solved.  Selecting  the  Solve  button  will  start  the  execution  of  the 
optimization  code.  The  SCTB  uses  MINOS  5.4  to  solve  the  optimization  problem  [14]. 
Depending  on  the  number  and  type  of  the  constraints  and  the  number  and  type  of  the 
basis,  the  optimization  can  be  solved  very  quickly  or  very  slowly.  Through  experience 
with  solving  SISO  problems  for  systems  of  8th  order  or  less,  it  is  determined  that  if  the 
basis  size  is  under  100,  the  constraints  are  feasible  for  a  basis  within  100  functions, 
and  the  SCTB  is  run  on  a  Sparc  20  system  or  better,  the  user  can  expect  typical 
MINOS  runs  to  take  less  than  ten  minutes. 

When  MINOS  starts  solving  the  problem,  it  starts  producing  output  to  the  MAT- 
LAB  workspace.  There  are  two  steps  to  a  the  MINOS  solution  process.  At  first, 
MINOS  tries  to  find  a  feasible  solution.  Here  the  columns  of  interest  in  the  MINOS 
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Figure  5.13:  Constrained  Optimization  Control  Design  Panel:  Post-Solution 

output  are  the  Itn,  ninf ,  and  sinf  columns.  The  Itn  column  is  the  iteration  num¬ 
ber.  For  a  paxticulax  iteration,  the  ninf  column  shows  the  number  of  infeasibilities 
and  the  sinf  column  shows  the  sum  of  those  infeasibilities.  As  MINOS  reaches  a 
feasible  solution,  the  ninf  column  and  the  sinf  should  decrease  to  zero,  although 
the  sinf  column  may  have  slight  increases. 

Once  the  feasible  solution  is  reached,  the  second  portion  of  the  solution  process 
aims  to  optimize  the  objective  function.  Here,  the  objective  column  becomes  a 
column  of  interest.  This  is  the  objective  value,  which  MINOS  is  working  to  minimize, 
for  a  particular  iteration. 

It  is  possible  that  there  is  no  feasible  solution.  When  this  occurs,  the  number 
and  sum  of  the  infeasibilities  may  show  how  close  the  problem  is  to  feasibility.  If 
the  sum  of  infeasibilities  is  very  small,  an  addition  of  more  basis  functions  or  a  slight 
relaxation  of  constraints  may  make  the  problem  feasible. 

If  the  problem  is  feasible,  then  a  controller  will  be  sent  to  the  SCTB  and  the 
Constrained  Optimization  Controller  Design  panel  will  look  like  Figure  5.13.  There, 
the  order  of  the  controller  design  and  the  length  of  time  the  SCTB  used  to  solve  the 
optimization  problem  are  displayed.  The  final  objective  function  value  is  found  in  the 
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Figure  5.14:  Analysis  Tool:  Initial  Display 


MINOS  workspace  output. 

After  the  feasible  solution  is  found,  the  controller  needs  to  be  analyzed  to  de¬ 
termine  if  it  is  a  satisfactory  solution.  A  controller  that  meets  all  constraints  may 
not  be  a  desirable  controller.  There  may  be  a  spiking  problem,  where  the  solution 
leaks  through  holes  in  constraints.  Spiking  is  corrected  by  increasing  the  density  of 
the  constraints.  Settling  problems  might  also  occur  in  the  time  response,  settling 
constraints  are  not  extended  to  infinity.  Depending  on  when  and  where  the  settling 
constraints  end,  the  solution  may  not  have  actually  settled.  In  this  case,  the  solution 
may  be  unstable,  or  may  stabilize  at  a  time  much  later  that  desired.  The  remedy  is 
to  extend  the  settling  constraints  to  a  further  point  in  time. 

5.4  Analysis  Tool 

To  invoke  the  Analysis  Tool,  select  the  Time/FVeq  Analysis  item  of  the  Analysis 
menu  bax.  The  initial  display  of  the  Analysis  Tool  is  shown  in  Figure  5.14.  In  the 
Analysis  Tool,  there  are  four  buttons  in  the  Analysis  Options  section:  1  Plant 
w/Multiple  Controllers,  1  Controller  w/Multiple  Plants,  Multiple  Plants, 
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Figure  5.15:  Analysis  Tool:  Closed  Loop  Display  Option 

and  Multiple  Controllers.  The  1  Plant  w/Multiple  Controllers  button  is 
depressed  in  the  initial  display  with  the  corresponding  menus  for  selecting  a  single 
plant  and  multiple  controllers.  There  are  two  buttons  in  the  Display  Options  section: 
Loop  Transfer  Function  and  Closed  Loop.  The  Loop  Tiransfer  Function 
button  is  selected  in  the  initial  display.  There  axe  also  four  other  buttons  in  the 
Analysis  Tool  GUI:  Refresh,  LTI  View,  View,  and  Close. 


5.4.1  Analysis  Options 

The  first  of  the  Analysis  Options  is  the  1  Plant  w/Multiple  Controllers  option. 
If  the  Loop  IVansfer  Function  button  is  selected,  the  Analysis  Tool  will  prepare  to 
plot  open  loop  transfer  functions  of  the  plant /controller  combinations  of  the  selected 
plant  and  controllers.  The  selected  controllers  must  have  the  same  number  of  inputs 
and  outputs.  If  the  Closed  Loop  button  is  selected,  the  input /output  selection  grid 
and  list  of  weighting  functions  are  displayed  as  in  Figure  5.15.  The  input/output  pair 
usually  selected  is  the  djc  pair.  For  the  SISO  loops  in  this  thesis,  the  (<i,c)  pair  will 
contain  all  the  outputs  (defined  as  e  outputs)  in  relation  to  the  commanded  input 
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(defined  as  a  d  input).  The  weighting  options  are  not  used  in  this  thesis. 

The  second  of  the  Analysis  Options  is  the  1  Controller  w/Miiltiple  Plants 
option.  This  option  works  just  like  the  1  Plant  w/Multiple  Controllers  option. 
The  only  difference  is  that  there  are  multiple  plants  selected  for  one  controller.  The 
selected  plants  all  need  to  have  the  same  inputs  and  outputs.  The  third  of  the 
Analysis  Options  is  the  Multiple  Plants  option,  for  comparing  plants  to  other 
plants.  When  selected  by  itself,  the  Multiple  Plants  option  will  automatically  have 
the  input /output  selection  grid  and  list  of  weighting  functions  displayed  while  the 
Display  Options  section  disappears.  The  Multiple  Plants  option  may  be  used  in 
conjunction  with  the  1  Plant  w/Multiple  Controllers  option  or  the  1  Controller 
w/Multiple  Plants  option.  In  either  case,  the  Display  Options  section  reappears. 
The  final  option  is  the  Multiple  Controllers  option,  which  compares  controllers  to 
other  controllers.  In  this  case,  there  are  no  Display  Options,  no  input /output  grid, 
and  no  weighting  function  options. 

5.4.2  Analysis  Plots 

Since  all  the  plants  and  controllers  will  have  the  same  sample  rate  (for  discrete  sys¬ 
tems)  and  are  not  FRF  objects,  the  LTI  View  option  is  the  only  option  necessary 
(making  View  an  unnecessary  option).  LTI  View  will  combine  the  selected  systems 
and  call  the  MATLAB  Itiview  function.  This  will  invoke  the  LTI  Viewer.  An  im¬ 
portant  feature  of  the  LTI  Viewer  is  the  Plot  Type  menu  bar.  There  are  many  plot 
options,  but  the  two  plot  options  used  will  be  the  Bode  and  Step,  which  will  plot 
the  Bode  Plots  and  Step  Responses  respectfully.  The  Zoom  bar  allows  for  horizontal, 
verticcJ,  or  combination  zooming.  There  is  2ilso  a  Full  View  button  that  returns  the 
zoomed  plot  to  its  original  dimensions. 

There  are  foiur  pull-down  menus:  File,  Tools,  Plots,  and  Help.  The  File  menu 
has  three  options:  New  Viewer,  Print  Response,  and  Close  Viewer.  New 
Viewer  creates  a  new  viewer  window  and  Close  Viewer  closes  the  viewer.  Print 
Response  allows  the  user  to  print  the  plot  on  the  screen  to  either  a  printer  or  to  a  file. 
Figure  5.16  is  an  example  of  printed  LTI  Viewer  output.  The  Tools  menu  has  three 
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Figure  5.16:  LTI  Viewer  Output 

options:  Viewer  Controls,  Response  Preferences,  and  Linestyle  Preferences. 
Viewer  Controls  toggles  between  having  viewer  controls  available  or  displaying  the 
plot  by  itself.  Response  Preferences  controls  the  axes  and  units  of  the  plots. 
Linestyle  Preferences  controls  the  lines  of  the  plots.  Different  linestyles  include 
different  colored  lines,  dashed  and  dotted  lines,  and  lines  conaposed  of  repeated  sym¬ 
bols.  The  Plots  menu  has  a  solitary  Grid  On  option  that  toggles  between  displaying 
gridlines  on  the  plots.  Help  is  standard  and  self  explanatory. 


5.5  Model  Reduction 

To  reduce  a  model,  select  the  controller  (or  plant)  to  be  reduced,  and  select  the 
Reduction  item  of  the  respective  Pre-Processor  Main  Panel  menu  bar.  This  will 
invoke  the  Model  Reduction  Tool  as  in  Figure  5.17.  To  reduce  the  model,  a  method 
must  first  be  selected  out  of  the  Method  menu  bar.  There  are  many  different  meth¬ 
ods  available,  but  this  thesis  concentrated  on  the  Balance  and  Truncate  and  the 
Fractional  Balanced  Reduction  methods,  both  of  which  are  offered  in  the  Method 
menu  bar. 
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Figure  5.17:  Model  Reduction  Tool 


Once  the  method  is  determined,  selecting  the  Reduce  button  will  invoke  a  fre¬ 
quency  weighting  option.  When  this  option  is  accepted,  the  Frequency  Weighting 
Tool  is  invoked.  There  axe  options  for  a  lowpass,  highpaiss,  or  notch  filter.  This  filter 
can  be  placed  at  the  input,  output,  or  both.  After  the  filter  is  done  (or  if  there  was 
no  weighting),  the  Reduction  Tool  is  invoked  as  in  Figure  5.18.  This  tool  allows  the 
user  to  select  how  many  modes  the  reduced  order  model  will  have.  The  plot  on  the 
Reduction  Tool  provides  a  criterion  for  selecting  modes.  Note  that  in  Figure  5.18, 
the  minimum  normalized  Hankel  singular  value  over  the  fifth  most  significant  mode 
is  similar  to  that  for  the  25th  most  significant  mode.  This  means  that  the  five  mode 
reduced  system  is  almost  as  good  the  25  mode  reduced  system. 

Once  the  model  is  reduced,  the  reduced  order  model  is  placed  in  the  System  menu 
of  the  Model  Reduction  Tool.  This  system  can  be  saved  as  a  controller  in  the  Main 
Panel  by  hitting  the  Save  button  in  the  Model  Reduction  Tool. 


5.6  The  Close  Loop  Tool 

The  final  SCTB  application  to  be  discussed  is  the  Close  Loop  Tool.  Before  the  Close 
Loop  Tool  can  be  used,  the  active  plant  and  active  controller  must  be  selected.  To 
invoke  the  Close  Loop  Tool,  shown  in  Figure  5.19,  select  the  Close  Loop  button 
in  the  Main  Panel.  Select  the  Close  Loop  button  in  the  Close  Loop  Tool  GUI  to 
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Figure  5.18:  Reduction  Tool 


Figure  5.19:  Closed  Loop  Tool 
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Figure  5.20:  Closed  Loop  Tool:  Save  Panel 

close  the  loop  around  the  active  controller.  This  will  invoke  the  Save  Panel.  This 
panel,  shown  in  Figure  5.20  can  show  the  Bode  plots,  if  desired.  This  window  will 
also  always  state  that  the  closed  loop  system  is  unstable.  Closed  loop  stability  should 
be  determined  before  closing  the  loop.  To  save  the  closed  loop  plant,  select  the  Save 
button  in  the  Save  Panel.  This  will  save  the  closed  loop  system  as  a  plant  in  the 
Main  Panel  and  close  all  associated  Close  Loop  Tool  windows. 
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Chapter  6 


Draper  Small  Autonomous  Aerial 
Vehicle  Control  Example 


The  Draper  Small  Autonomous  Aerial  Vehicle  (DSAAV)  is  a  small  model  helicopter 
equipped  with  an  avionics  system  designed  by  a  team  from  MIT,  Boston  University, 
and  the  Draper  Laboratory  for  the  1996  International  Aerial  Robotics  Competition. 
Using  classical  control  techniques,  the  DSAAV  successfully  performed  completely  au¬ 
tonomous  missions  and  won  the  competition.  Since  the  contest,  the  DSAAV  has 
become  a  research  platform  for  advanced  flight  control,  among  other  things.  An 
example  of  this  is  described  in  [16]. 

6.1  Design  Objectives 

The  primary  design  objective  is  to  improve  the  performance  of  the  design  through 
reduction  of  the  error  responses  of  the  inner  loops.  The  classical  controllers  previ¬ 
ously  designed  for  the  DSAAV  are  the  standards  to  whidi  this  design  will  be  held. 
Another  major  objective  is  to  determine  the  utility  of  decoupling  the  channels  into 
the  hierarchical  structure,  as  opposed  to  leaving  each  channel  coupled  as  a  MIMO 
system.  A  further  design  objective  is  to  determine  the  effect  of  the  constraints  on  the 
solutions. 

There  are  also  objectives  within  the  procedure.  One  is  to  determine  the  tradeoff 
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Figure  6.1:  Entire  DSAAV  System  of  Coupled  Channels 

between  using  the  Basis  Function  Algorithm  and  finding  a  high  order  solution  through 
a  large,  inefficient  basis.  It  is  also  of  interest  to  determine  how  well  the  time  scale 
iteration  of  the  Basis  Function  Algorithm  works.  A  final  objective  is  to  determine 
the  utility  of  the  Q-minimization  objective  function. 

6.2  Decoupled  Hierarchical  System  Architecture 

The  main  DSAAV  system  is  depicted  in  Figure  6.1.  There  are  four  main  channels  for 
each  direction  of  motion:  forward  position  r,  lateral  position  y,  vertical  position 
and  heading  ip.  Each  channel  has  an  associated  primary  control  actuator.  Forward 
motion  has  the  main  rotor  pitch  cyclic  i7(9cyc?  Istteral  motion  has  the  main  rotor  roll 
cyclic  C/^cyc,  vertical  motion  has  the  main  rotor  collective  UmrcoI^  and  heading  has 
the  tail  rotor  collective  UtrcoI*  Since  each  channel  has  an  associated  control  actuator, 
it  is  assumed  that  each  channel  can  have  its  own  associated  controller. 

The  dynamics  of  the  four  channels  are  coupled.  The  first  decoupling  separates 
the  dynamics  of  the  entire  system  into  the  four  channels.  More  coupled  dynamics 
exist  within  each  channel.  The  coupled  forward  motion  MIMO  channel  is  depicted  in 
Figure  6.2.  The  pertinent  states  for  the  forward  motion  channel  are  forward  position 
re,  forward  velocity  w,  pitch  6^  and  pitch  rate  q.  Each  channel  is  decoupled  into 
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Figure  6.2:  Coupled  MIMO  Structure  for  DSAAV  Forward  Motion 

the  hierarchical  architecture  of  SISO  loops,  as  in  Figure  6.3  for  the  forward  motion 
channel.  Here  each  pertinent  state  is  given  its  own  SISO  loop,  complete  with  its  own 
commanded  state,  dynamics,  and  controller.  Each  SISO  loop  has  its  own  commanded 
state  as  the  primary  loop  input.  This  commanded  state  is  the  output  of  the  previous 
loop’s  controller.  For  example,  the  commanded  pitch  rate  qc  is  the  output  of  the  pitch 
controller  Ke. 

6.2.1  Forward  Motion  Dynamics 

The  mathematical  decoupling  of  the  forward  motion  starts  with  the  state  space  de¬ 
scription  of  the  coupled  dynamics  Pfwd{s)-  There  are  two  internal  states  for  the 
forward  motion  loop:  the  main  rotor  plane  pitch  angle  q,\  and  the  flybar  plane  pitch 
angle  ci/j,.  The  state  space  representation  of  the  coupled  system  Pfwd{s)  is: 

x  =  Ax  +  BUecyc  (6.1) 
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with  the  states  being  x  =  [xu6qai 


(6.2) 


6.2.2  The  Pitch  Rate  Loop 

The  innermost  loop  of  the  forward  motion  channel  is  the  pitch  rate  loop.  To  construct 
the  pitch  rate  loop,  the  pitch  rate  plant  P^is)  is  created.  P4{s)  is  a  truncated  version 
of  Pfwdi  with  the  dynamics  associated  with  x,  u,  and  $  completely  removed.  The 
pitch  rate  loop  is  depicted  in  Figure  6.4.  The  input  is  the  commanded  pitch  rate  qc. 
The  resulting  pitch  rate  plant  P4{s)  is  the  third  order  system: 

7154s +  91271  _  7154(s  +  12.82) _ 

“  s3  +  53.903s2  +  323s  +  5285  (s  +  2.1834  ±  10.096i)(s  +  49.536) '  ^  ^ 

The  state  space  representation  of  the  model  based  pitch  rate  plant  Pq{s),  as  de¬ 
scribed  in  Figure  6.5,  has  the  truncated  states  =  [q  di  This  state  space 

system  is  as  follows: 

Xq  —  i4g®^  +  BqV^ 

Uq  =  Cq^  +  DqUq 
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(6.4) 

(6.5) 


Figure  6.5:  Pitch  Rate  Loop:  Model  Based  Representation 


with  the  state  space  matrices  defined  as: 


■^g  ®oi)9  ®oi,oi  ®oitOl/6  ^ 

^aifb,g  ^aifbiOi  ^aifb,aifb  ^oi/fc  ® 

-100  01 

(7,=  100  -D,  =  00 

0  0  0  1  0 


(6.6) 


with  the  inputs  =  [Uecyc  QcV  being  the  control  aetuator  Uecyc  and  the  commanded 
pitch  rate  qc.  The  outputs  ^  =  [?err  q  Uecyc]'^  are  the  pitch  rate  error  gerr  (which 
is  fed  into  controller  Kq),  the  pitch  rate  itself,  and  the  control.  These  outputs  were 
necessary  for  analysis  in  the  SCTB. 


Closing  the  Loop  Following  Pitch  Rate  Controller  Design 

Once  Pq{s)  is  defined,  the  pitch  rate  controller  Kq  is  ready  for  design  using  the 
solution  procedure.  This  design  is  described  later  in  this  chapter.  The  state  space 
representation  of  Kq  is  as  follows: 


Xj^q  —  Ajcq^^q  "t"  B^q(Jerr 

BkqQ  "1"  Bkq^c  ”1"  “^kq^q 

(6.7) 

W  “  ^kq^^q  "i“  ^kqQerr 

~^kqQ  J^kqQc  H“  ^kq^^q 

(6.8) 

where  xj^q  are  the  controller  states  of  Kq  and  Akg  is  a  Ukq  x  Ukq  matrix. 
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After  Kq  is  designed,  close  the  loop  to  obtain  the  closed  loop  system  Hqi 

^  +  BqQc 


(6.9) 


^  =  Cqg^  +  Dqqc 

with  the  state  space  matrices  defined  as: 


Aq  = 


Aq{:,l)-Bq{:,l)Dkq  A,(:,2:3)  Bq{:,l)Ckq 
~Bkq  0ix2  Akq 

1  0  0  0 
Dkq  0  0  Ckq 


Cq  = 


(6.10) 


Dq  = 


Bkq 

0 

Dkq 


(6.11) 

where  the  states  are  ^  ^q]^  and  the  outputs  are  =  [?  UocycV •  The  input 

is  the  commanded  pitch  rate. 


6.2.3  The  Outer  Three  Loops 

Construction  of  the  pitch,  forward  velocity,  and  forward  position  loops  is  very  similar 
to  the  pitch  rate  loop  construction.  The  construction  of  these  loops  used  the  proce¬ 
dures  described  in  Chapter  3  for  intermediate  and  outer  loops.  Therefore,  this  section 
will  not  detail  how  the  outer  three  loops  were  constructed. 


6.3  Controller  Design  using  Solution  Procedure 

The  standard  to  which  the  controller  designs  will  be  measured  against  is  the  classical 
controller  designs  presently  used  on  the  DSAAV.  The  DSAAV  is  a  discrete  time 
system  with  a  sample  rate  of  25  Hertz  (or  0.04  seconds).  The  current  actual  classical 
designs  for  pitch  rate  KqA  (the  ’A’  in  KqA  stands  for  actual)  forward  velocity  KuA  are 
first  order  controllers.  The  current  actual  classical  designs  for  pitch  Koa  and  forward 
position  KxA  are  proportional  controllers. 

The  performance  standard  to  which  the  controller  designs  will  be  held  in  this 
section  are  the  closed  loop  step  responses  (based  on  a  unit  step  of  the  commanded 


95 


step  Resxxise 

'^FromT*  (j. 


Figure  6.6:  Pitch  Rate  Error  Step  Response:  Based  on  Controller  KqA 

state)  of  the  state  errors.  For  the  pitch  rate  loop,  this  would  be  the  pitch  rate  error 
transient  qerr{t)-  The  step  responses  for  controllers  KqA,  K$a,  KuA,  and  Kxa  are  in 
Figures  6.6  through  6.9  respectively. 

The  rest  of  this  section  will  detail  the  application  of  the  solution  procedure.  The 
pitch  rate  loop  controller  {Kq)  design  will  be  used  as  the  detailed  example  for  the 
solution  procedure.  Only  the  resulting  solutions  to  the  other  three  loops  will  be 
presented. 

6.3.1  Pitch  Rate  Loop:  Constrained  Optimization  Prepara¬ 
tion 

This  subsection  sets  up  the  constrained  optimization  problem  for  the  pitch  rate  loop 
controller.  This  preparation  describes  the  final  objective  function  and  constraints 
used  to  find  the  best  solution. 

The  initial  objective  function  is  always  the  'H2  minimization  objective.  This  ob¬ 
jective  function  was  also  used  to  find  the  pitch  rate  loop  controller.  For  this  entire 
example,  a  first  order  LPF  with  a  10  Hertz  cutoff  frequency  was  always  added  to  the 
input.  Notice  that  10  Hertz  is  less  than  half  of  the  system  sample  rate. 

Throughout  the  solution  procedure,  the  constraints  selected  were  based  on  the 
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Figure  6.7:  Pitch  Error  Step  Response:  Based  on  Controller  Kba 
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Figure  6.8:  Forward  Velocity  Error  Step  Response:  Based  on  Controller  KuA 


StepR^™ 


Figure  6.9:  Forward  Position  Error  Step  Response:  Based  on  Controller  Kxa 

closed  loop  characteristics  of  KqA^  The  constraints  presented  were  those  used  to 
find  the  design  solution  and  were  altered  from  the  initial  constraints  as  part  of  the 
iterative  nature  of  the  solution  procedure.  The  frequency  constraint  for  the  sensitivity 
transfer  function  iS'(5)  (for  command  following)  are  the  points  (0.00001  Hz,  —50  dB), 
(0.001  Hz,  -40  dB),  (0.01  Hz,  -20  dB),  (0.1  Hz,  -5  dB),  and  (2  Hz,  1.25  dB).  The  first 
four  point  constraints  were  chosen  to  push  the  low  to  middle  frequency  gain  further 
than  KqA,  for  command  following.  The  last  constraint  was  chosen  to  control  the  high 
frequency  peaJc.  The  sensitivity  transfer  function  constraints  with  the  corresponding 
closed  loop  curve  for  controller  KqA  are  shown  in  Figure  6.10. 

The  frequency  constraint  for  the  closed  loop  transfer  function  C7(s)  (for  stability 
robustness)  consists  of  two  line  constraints.  One  is  an  upper  bound  line  at  0  dB 
from  0.01  Hertz  to  2  Hertz.  The  other  is  the  roll-off  line  from  (1  Hz,  0  dB)  to 
(10  Hz,  —40  dB).  These  line  constraints  for  the  closed  loop  transfer  function  and  the 
corresponding  closed  loop  curve  for  controller  KqA  are  shown  in  Figure  6.11. 

After  much  experience  with  the  SCTB  Constrained  Optimization  Tool,  it  is  deter¬ 
mined  that  a  settling  constraint  is  often  sufficient  for  a  time  constraint.  If  the  settling 
constraint  is  not  too  aggressive,  the  amplitude  of  the  time  response  will  remain  within 
limits.  Aggressive  settling  constraints  will  result  in  solutions  that  compensate  with 
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Figure  6.10:  S{s)  Frequency  Constraint  Based  on  KqA 


Input  (^;  Output  q 


Figure  6.11:  C'(s)  Frequency  Constraint  Based  on  KqA 
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Figure  6.12:  qerr{t)  Step  Response  Constraint 


very  high  error  amplitude  peaks  before  meeting  the  settling  constraint.  The  time 
constraint  for  the  step  response  5err(t)  is  to  settle  within  il%  error  from  8  to  15  sec¬ 
onds.  This  settling  constraint  and  the  corresponding  step  response  for  controller  KqA 
is  shown  in  Figure  6.12. 


6.3.2  Pitch  Rate  Loop:  Basis  Function  Algorithm 

Once  the  objective  and  constraints  axe  in  place,  the  Basis  Function  Algorithm  can  be 
applied  to  find  an  efficient  solution.  Based  on  experience  and  the  size  of  the  problem, 
the  toleration  limit  on  the  number  of  basis  functions  Ntoi  is  set  to  100.  Also,  based 
on  the  size  of  the  problem,  the  acceptable  bound  separation  range  A  is  set  to  1. 

Two  cases  of  the  Basis  Function  Algorithm  were  applied.  The  first  case  did  not  use 
the  time  scale  iteration  for  Legendre  and  Laguerre  functions.  Instead,  the  time  scale 
frequencies  of  1.25  Hertz  for  Legendre  functions  and  6.25  Hertz  for  Laguerre  functions 
were  applied.  These  were  chosen  because  they  are  the  nominal  frequencies  that  the 
SCTB  will  suggest  if  a  specific  time  scale  is  not  entered.  This  Laguerre  frequency  is 
one  fourth  of  the  sampling  rate  and  this  Legendre  frequency  is  one  twentieth  of  the 


sample  rate.  This  is  the  only  reasoning,  other  than  convenience,  behind  specifically 
choosing  these  frequencies  over  any  other.  The  other  case  used  the  time  scale  iteration 
to  try  to  come  up  with  an  efficient  time  scale  for  Legendre  and  Laguerre  functions. 

As  a  result  of  numerous  iterations  (involving  altering  constraints)  of  the  Basis 
Function  Algorithm,  an  acceptable  solution  was  found  after  the  final  set  of  constraints 
was  determined  to  be  those  in  the  previous  subsection.  This  best  designed  controller 
KgD32  (the  ’D’  stands  for  design),  was  considered  acceptable  when  its  first  order 
reduction  passed  implementation  on  the  DSAAV  simulator.  KqD32  was  found  using 
the  fixed  time  scale  case.  The  corresponding  basis  was  comprised  of  25  Laguerre 
functions  at  6.25  Hertz,  two  Legendre  functions  at  1.25  Hertz,  and  two  FIR  functions. 
This  produced  a  solution  controller  of  32nd  order.  To  reach  this  solution,  the  Basis 
Function  Algorithm  solved  40  optimization  problems,  each  between  45  seconds  and 
75  seconds. 

When  the  Basis  Function  Algorithm  with  time  scale  iterations  was  applied  to  the 
same  constraints,  it  produced  an  efficient  8th  order  solution  controller.  The  basis  was 
comprised  of  four  Laguerre  functions  at  0.5  Hertz  and  one  FIR  function.  To  reach 
this  solution,  the  Basis  Function  Algorithm  solved  71  optimization  problems,  each 
requiring  between  17  seconds  and  64  seconds.  The  first  order  reduction  appeared  to 
be  acceptable  in  the  SCTB  analysis  phase.  However,  it  had  large  pitch  oscillations  in 
the  DSAAV  simulator  test.  Therefore,  this  was  an  unacceptable  design. 

6.3.3  Pitch  Rate  Loop:  Model  Reduction  and  An^ysis 

After  DSAAV  simulator  test  failures  for  all  controllers  of  third  order  and  higher,  it  was 
determined  that  controllers  that  are  second  order  or  less  are  implementable.  This  led 
to  the  model  reduction  of  R',r>32.  The  FBR  model  reduction  method  reduced  KqD32 
into  the  first  order  controller  Kqo- 


f  .\  0.02481 («+2.3729)  .  \  0.02598^-0.02363 

—  «+0.002497  ^qD\^)  z-0.9999 


(6.12) 


It  is  noted  that  this  design  is  very  similar  to  a  proportional  integral  (PI)  controller. 
The  sensitivity  transfer  function  curves  for  KqA  and  Kqo  are  shown  in  Figure  6.13. 
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Figure  6.13:  S{s)  for  Pitch  Rate  Controller  Designs 


Here,  Kqo  (the  darker  curve)  starts  at  -50  dB  for  the  low  frequency  range,  suggesting 
good  command  following.  At  the  high  frequency  end,  the  two  designs  share  similar 
dynamics. 

The  closed  loop  transfer  function  curves  for  the  controller  designs  are  shown  in 
Figure  6.14.  Here,  both  designs  roll  off  well  at  the  high  frequency  end.  Although  the 
classical  controller  (darker  curve)  is  at  a  lower  gain  than  the  constrained  optimization 
controller,  both  curves  roll  off  well  enough  to  suggest  sufficient  stability  robustness. 

The  step  responses  for  the  controller  designs  are  shown  in  Figure  6.15.  The  step 
response  for  the  constrained  optimization  controller  (solid  line)  performed  well,  set¬ 
tling  to  ±1%  error  in  6.4  seconds.  This  represents  a  significant  improvement  over  the 
classical  design,  which  settles  to  ±1%  error  in  81.1  seconds. 


DSAAV  Simulator 

To  determine  the  implementabihty  of  the  designs,  they  were  tested  on  a  DSAAV 
simulator.  This  simulator  shows  a  picture  of  the  DSAAV  helicopter  in  a  field  as 
it  flies  the  predetermined  flight  script  on  the  screen.  There  are  four  main  possible 
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Figure  6.15:  Step  Responses  for  Pitch  Rate  Controller  Designs 
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results  from  a  simulator  flight  test.  One  is  a  smooth  flight  where  the  DSAAV  flies  its 
script  well  without  oscillations  or  diversions  from  its  flight  path.  The  second  result 
is  that  the  DSAAV  will  oscillate  throughout  its  entire  flight.  For  the  forward  motion 
loop,  these  will  be  pitch  oscillations.  If  the  oscillations  are  moderate  enough,  the 
DSAAV  will  fly  and  complete  its  flight  script,  oscillating  most  of  the  way.  Oscillations 
represent  a  stability  problem.  The  third  result  is  the  DSAAV  will  deviate  from  its 
scripted  flight  path.  In  the  moderate  case,  the  control  does  not  follow  commands 
well  and  passes  the  scripted  landing  point.  A  higher  level  guidance  loop  will  try 
to  compensate  for  this  and  turn  the  DSAAV  around  in  an  attempt  to  return  to 
the  scripted  flight  path.  The  extreme  case  has  the  DSAAV  immediately  deviating 
from  the  entire  flight  path.  This  is  usually  coupled  with  oscillations  and  represents  a 
stability  and  command  following  problem.  The  final  case  is  where  the  DSAAV  crashes. 
In  the  simulator,  when  the  DSAAV  crashes,  it  will  often  bounce  off  the  ground  and 
flip  around  like  a  football  until  the  simulator  decides  to  end  the  simulation  with  an 
error  message.  The  crashing  DSAAV  indicates  a  design  that  is  poor  and  can  have 
problems  that  may  not  be  apparent.  At  the  very  least,  it  is  obvious  that  there  is  a 
stability  problem. 

6.3.4  Solution  Controllers  of  Outer  Loops 

After  the  design  of  Kqo,  the  solution  procedure  was  applied  to  the  pitch  loop.  The 
constrained  optimization  tool  was  able  to  design  higher  order  controllers  that  were 
improvements  upon  K$A‘  Unfortunately,  when  these  controllers  were  reduced  to  an 
implementable  order,  the  performance  loss  made  the  reduced  order  designs  worse 
than  Koa-  Since  the  classical  controller  was  just  a  proportional  controller,  it  was 
decided  to  keep  Kqa  as  the  pitch  controller.  There  is  still  a  difference  between  the 
classical  design  and  the  constrained  optimization  design  from  the  difference  in  pitch 
rate  controllers.  This  is  depicted  in  the  closed  loop  step  response  of  the  pitch  error  in 
Figure  6.16.  The  constrained  optimization  design  (solid  line)  has  more  overshoot  than 
the  classical  design,  but  reduces  the  settling  time  from  13.7  seconds  to  4.2  seconds. 

The  forward  velocity  design  was  very  similar  to  the  pitch  rate  design.  The  reduced 
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Figure  6.16:  6err{t)  for  Pitch  Controller  Designs 
first  order  design  K^d  is  as  follows; 

T,  /.N  -0.06531(*+0.078956)  r>-  /.A  -0.06541z+0.06521 

—  a+0.001837  ^uD\^)  *-0.9999  ’  '  ' 

This  was  reduced  from  a  37th  order  design.  It  is  also  noted  that  KuD  is  similar  to 
a  PI  controller.  The  comparison  between  the  classical  design  and  the  constrained 
optimization  design  is  shown  in  the  closed  loop  step  response  of  the  forward  velocity 
error  in  Figure  6.17.  The  constrained  optimization  design  (solid  line)  smoothes  the 
error  curve  reducing  the  classical  design  overshoot  from  41%  to  30%  while  maintaining 
a  similar  settling  time  (11.6  seconds  to  13.8  seconds).  Most  importantly,  the  forward 
velocity  design  passed  the  DSAAV  simulation  test. 

After  the  design  of  it  was  time  to  design  the  outermost  forward  position  loop 
controller  Kxd-  Since  the  classical  controller  Kxa  was  a  proportional  controller,  like 
KeA,  it  was  decided  to  keep  Kg  a  as  the  forward  position  controller.  The  result  was  a 
miniTnal  difference  between  classical  design  and  the  constrained  optimization  design. 
This  is  shown  in  the  closed  loop  step  response  of  the  forward  position  in  Figure  6.18. 
The  only  difference  between  the  designs  is  not  very  remarkable.  This  difference  is 
that  the  constrained  optimization  design  (solid  line)  reduced  the  settling  time  from 
51.9  seconds  to  51.8  seconds. 
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Figure  6.17:  tterr(i)  for  Forward  Velocity  Controller  Designs 
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Figure  6.18:  Xerrit)  for  Forward  Position  Controller  Designs 
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6.3.5  Controller  Design  for  Coupled  Forward  Motion 

To  compare  the  controller  design  process  of  the  SISO  hierarchical  structure  to  the 
entire  coupled  MIMO  system,  the  same  constraints  that  were  applied  to  the  SISO 
solutions  in  the  previous  subsections  were  applied  to  the  coupled  MIMO  forward 
motion  channel.  Not  only  was  the  SCTB  unable  to  find  a  feasible  solution,  it  was 
taking  well  over  two  hours  to  determine  that  the  problem  was  infeasible.  In  addition, 
the  SCTB  was  returning  unsolved  controllers  of  over  600th  order.  If  there  was  a 
feasible  solution  to  the  coupled  problem,  it  would  take  way  too  long  and  be  way 
too  complex  to  be  beneficial.  With  the  hierarchical  structure,  the  solutions  were 
always  found  in  less  than  ten  minutes  and  were  rarely  solved  in  more  than  two  and 
a  half  minutes.  This  allowed  the  Basis  Function  Algorithm  to  solve  the  optimization 
problem  40  to  80  times  for  one  design. 


6.3.6  (^-Minimization  Controller  Design 

The  constrained  optimization  controllers  designed  for  the  SISO  loops  of  the  DSAAV 
forward  motion  channel  were  all  based  on  the  'H2  objective  function.  Although  the 
'H.2  objective  was  sufficient  in  designing  these  controllers,  the  Q-minimization  ob¬ 
jective  required  investigation.  Using  the  designed  pitch  rate  controller  KqD  as  the 
nominal  controller,  the  ^-minimization  objective  was  used  to  design  a  pitch  rate  con¬ 
troller  to  modestly  improve  upon  the  Kqo  design.  The  same  frequency  constraints 
used  to  design  Kqo  were  applied  to  the  (^-minimization  design.  The  difference  was  a 
slight  tightening  of  the  step  response  constraint:  to  settle  within  ±1%  error  from  6 
to  15  seconds  (instead  of  8  to  15  seconds). 

The  Q-minimization  pitch  rate  controller  was  very  similar  to  Kqo-  The  reduced 
first  order  design  KqQ  is  as  follows: 


f  .\  0.022458(«+2.618)  v  _  0.02363z-0.02128 

~  *+0.0024945  *-0.9999 


(6.14) 


This  was  reduced  from  a  6th  order  design.  KqQ  is  also  very  similar  to  a  PI  con¬ 
troller  design.  The  improvement  was  very  modest:  reducing  the  settling  time  from 
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6.4  seconds  to  6.2  seconds.  Most  importantly,  KqQ  passed  the  DSAAV  simulation 
test. 

6.3.7  Observations  on  Constraints 

A  major  element  of  the  design  procedure  is  adjusting  constraints  to  achieve  feasible 
solutions  and  desired  performance.  Knowing  where  to  place  constraints  is  not  an 
exact  science  and  is  dependent  on  the  designer’s  experience  and  intuition.  Prom  this 
DSAAV  example,  there  were  noted  trends  regarding  constraints. 

For  time  constraints,  it  was  usually  sufficient  to  apply  a  settling  constraint.  This 
is  niiore  desirable  than  adding  an  amplitude  constraint  because  more  constraints  adds 
complexity  to  the  optimization,  resulting  in  infeasible  problems  or  higher  order  solu¬ 
tions.  When  applying  only  settling  constraints,  do  not  place  them  too  early.  Con¬ 
straints  that  force  the  error  to  decay  too  rapidly  will  result  in  error  curves  with 
high  amplitude  peaks  before  settling.  Settling  constraints  also  need  to  be  applied  for 
seven  to  twenty  seconds,  at  the  very  least.  Settling  constraints  that  do  not  endure 
can  lead  to  error  curves  that  either  do  not  settle,  settle  at  a  non-zero  value,  or  regain 
a  significant  amount  of  amplitude  before  settling. 

Often,  designs  that  performed  well  in  the  design  phase  did  not  perform  well  in  the 
DSAAV  simulator.  In  the  design  phase,  the  decoupled  dynamics  are  linearized,  where 
the  simulator  includes  the  non-linearized  dynamics.  Once  a  design  was  successful  in 
the  simulator,  the  frequency  domain  characteristics  of  that  design  can  be  used  for 
determining  the  implementability  of  future  designs.  For  a  given  loop,  implementable 
designs  (designs  that  pass  the  simulator)  often  would  have  similar  characteristics  in 
the  mid  to  high  frequency  range. 


6.4  Conclusions 

Replacing  the  pitch  rate  and  forward  velocity  classical  controllers  with  the  constrained 
optimization  designs  of  Kqu  and  KuD  showed  considerable  improvement  in  the  pitch 
rate  performance,  noticeable  improvement  in  the  pitch  and  forvvard  velocity  perfor- 
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mance,  and  negligible  improvement  in  the  forward  position  performance.  This  indi¬ 
cates  that  the  constrained  optimization  method  has  the  potential  to  find  controllers 
that  are  improvements  over  classical  designs.  In  addition,  the  Q-minimization  ob¬ 
jective  may  be  used  successfully,  in  the  procedure,  to  modestly  improve  upon  well 
designed  controllers. 

The  solution  procedure  was  able  to  find  successful  controllers.  However,  the  time 
scale  iteration  for  Legendre  and  Laguerre  functions  is  suspect.  As  shown  in  the  pitch 
rate  loop,  the  nominal  Legendre  and  Laguerre  time  scales  proved  to  be  more  successful 
than  the  detailed  time  scale  iteration.  This  is  despite  the  time  scale  iteration  finding 
more  efficient  (lower  order)  controllers.  Two  possible  explanations  come  to  mind. 
One  explanation  is  that  the  model  used  for  design  is  not  as  accurate  as  the  DSAAV 
simulator  model.  The  other  explanation  is  that  choosing  the  correct  time  scale  for 
these  functions  is  still  an  evolving  art,  requiring  further  investigation. 

Another  conclusion  is  that  the  hierarchical  structure  is  required  for  constrained 
optimization  of  highly  coupled  systems.  The  highly  coupled  MIMO  system  is  too 
complex  for  all  the  states  of  interest  to  be  constrained  at  the  same  time.  Although  the 
set  of  constraints  should  rem8un  feasible,  the  coupled  problem  is  too  computationally 
complex  for  the  constrained  optimization  procedure  to  solve.  Also,  the  probable 
solution  will  be  of  an  unmanageable  order.  The  hierarchical  structure  allows  the 
design  to  focus  on  each  pertinent  state  and  mahe  adjustments  swifter  and  easier. 
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Chapter  7 


Conclusions 


A  constrcdned  convex  optimization  control  design  procedure  that  takes  advantage  of  a 
decoupled  hierarchical  system  structure  has  been  presented  in  this  thesis.  The  solution 
procedure  was  applied  to  an  autonomous  helicopter,  to  show  a  practical  application  of 
the  procedure.  The  experimental  results  did  not  show  significant  improvement  in  the 
outer  loop  (position)  performance.  However,  within  the  coupled  channel,  there  was 
improvement  in  the  performance  of  the  inner  loop  dynamics  (attitude,  body  rates). 
Constrained  optimization  allowed  specifications  to  be  chosen  directly  as  constraints 
and  the  hierarchical  structure  allowed  for  controller  specialization,  focusing  on  im¬ 
provement  of  one  pertinent  state  at  a  time.  The  hierarchical  structure  also  allowed 
for  a  modular  design,  where  two  of  the  controllers  were  found  using  the  constrained 
optimization  solution  procedure  and  the  other  two  were  the  previously  designed  clas¬ 
sical  controllers.  Modular  design  is  attractive  because  different  approaches  might  be 
better  for  different  levels  of  the  hierarchy,  depending  on  the  dynamics  for  each  level’s 
control  problem. 

The  decoupled  hierarchical  structure  also  emphasized  a  significant  limitation  on 
constrained  optimization:  available  computational  power.  The  highly  coupled  MIMO 
control  system  is  almost  impossible  to  solve  with  constrained  optimization.  However, 
the  hierarchical  structure  makes  the  problem  tractable  by  breaking  the  problem  down 
into  simpler  SISO  loops.  Computer  weakness  is  still  an  issue  for  the  SISO  loops, 
hindering  constrained  optimization  if  there  are  too  many  constraints  or  too  many 
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basis  functions. 

The  constrained  optimization  solution  procedure  proved  to  be  effective  for  DSAAV 
controller  design.  In  addition,  the  procedure  demonstrated  the  advantages  of  the 
SCTB  MATLAB  linear  controller  design  tool.  However,  the  procedure  is  fax  from 
perfect.  The  procedure  was  able  to  estimate  an  appropriate  time  scale  for  the  Legen¬ 
dre  and  Laguerre  basis  functions  that  made  for  a  more  efficient  (lower  order)  solution. 
However,  these  efficient  solutions  were  proven  unimplementable,  failing  in  the  simula¬ 
tor.  This  is  in  contrast  to  higher  order  solutions  made  implementable  through  model 
reduction.  The  procedure  was  unable  to  distinguish  between  unimplementable,  effi¬ 
cient  solutions  and  higher  order  solutions  made  implementable  after  model  reduction. 

7.1  Recommendations  for  Future  Work 

The  most  difficult  portion  of  solving  the  constrained  optimization  control  problem  is 
finding  the  most  appropriate  basis  for  estimating  Q.  This  thesis  showed  an  example 
where  the  most  efficient  basis  (lowest  order)  for  a  set  of  constraints  may  not  be  the 
most  appropriate  basis  for  a  given  control  problem.  Therefore,  a  higher  order  solution 
based  on  a  higher  order  basis  might  provide  the  best  solution.  There  are  numerous 
design  variables  that  come  with  choosing  a  basis.  An  obvious  area  of  uncertainty  is 
deciding  how  many  of  each  type  of  basis  function  to  include.  Then  there  is  the  issue 
of  the  proper  time  scale  for  the  Legendre  and  Laguerre  basis,  if  they  are  selected.  If 
there  is  a  previous  design,  poles  from  that  Q  may  be  added  to  the  basis  as  well.  If  any 
of  these  decisions  axe  errant,  the  basis  may  be  inefficient  or  make  the  optimization 
problem  infeasible.  Therefore,  if  constrained  optimization  is  to  be  used  in  future 
control  problems,  as  it  should,  more  work  needs  to  be  done  to  demystify  the  art  of 
finding  the  best  basis.  At  the  very  least,  reducing  the  uncertainty  of  one  of  the  design 
variables  of  basis  choice,  such  as  time  scale,  can  greatly  simplify  the  problem. 

In  addition  to  the  tmcertainty  in  choosing  design  variables  for  the  basis  functions, 
there  axe  many  design  variables  when  it  comes  to  selecting  objective  functions  and 
constraints.  The  SCTB  gives  three  options  for  objective  functions.  This  thesis  only 
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fully  investigated  the  H2  minimization  objective  and  lightly  touched  upon  the  Q- 
minimization  objective.  The  objective  was  not  even  addressed.  In  addition,  the 
impact  of  the  frequency  weighting  on  the  objective  was  not  fully  addressed.  As 
for  constraints,  gain  margin  and  phase  margin  constraints  may  be  beneficial  to  a 
constrained  optimization  design. 

Another  avenue  for  future  investigation  is  a  more  systematic  procedure  for  map¬ 
ping  constraints  within  the  hierarchical  structure.  This  involves  determining  con¬ 
straints  on  the  innermost  loops  to  a  given  constraint  for  the  outermost  loop.  The 
objectives  of  the  outer  loops  would  then  depend  on  the  capabilities  of  the  lower  lev¬ 
els.  The  hierarchical  structure  should  provide  insight  into  determining  appropriate 
objectives  and  constraints  for  all  the  loops.  This  thesis  provided  some  intuitive  ob¬ 
servations  about  the  relationships  between  frequency  domain  constraints  and  closed 
loop  performance,  but  a  more  systematic  approach  is  needed  to  facilitate  automation 
of  the  constrained  optimization  procedure. 

As  mentioned  earlier,  computational  power  is  an  obstacle  to  constrained  optimiza¬ 
tion  control.  To  take  advantage  of  the  increasing  computing  power  that  comes  with 
the  advances  of  technology,  it  is  important  that  the  software  available  for  solving 
constrained  optimization  problems  is  always  adapting  and  improving  upon  the  latest 
advances  in  operations  research.  For  example,  MINOS  5.4,  the  software  used  to  solve 
the  constrained  optimization  problems  in  the  SCTB  is  dated  December  1992  and 
its  user’s  guide  is  dated  February  1995  [14].  Software  that  incorporates  the  model 
reduction  as  part  of  the  optimization  would  be  extremely  useful  for  systems  with 
implementation  limits. 
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